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(3) BE VE+EVG=0. 


As it is well-known that if we have R2,.,.=0, the space admits a group 
of affine motions of maximum order n?+m and vice versa, we consider only 
the case Rd,,.%+<0, 23, and assume r=7?2, for we shall shcw later (§ 3) 
that if Riv.» <0, 7 can not be more than 7? [I. P. Egorov, 1). As there are 
just 2 equations in (3), 7 can take on l,------ ,n, and we may treat 7 asa 
covariant index, for the equations in (3) may be substituted by linear combi- 
nations of them. 


$2. Curvature affinor of an A, admitting a group of affine motions of 
order 712, 


As the equations (2) must be linear combinations of (3), we get 


(4) Re eyo aaa Vow ’ 
< 5) R edk— ° you R weedy ° SD 
wr < A ee VW Ks. 


We first treat (5) to find the form of Rive. 

The first thing to do is to eliminate AP ves and we first consider the 
case where there is some vector u* such that the determinant | V2.2*| does 
not vanish, for we can prove later (§3) that if this determinant vanishes for 
any vector u*, we have R3,,.—0, and this contradicts the assumption 7=72. 

As | V2.w*|3=0, we can find out an affinor V® such that Vi,u*VP=88. 
Then multiplying (5) by w*V® and contracting, we get 


Ae Va Rive RR UVE—Rrrnw UV 8S 
—R oy uve, 
and substituting these into (5), we get equations of the form 
(6) ao Wve Ouse RR 
[ee ]w— ‘Re ue Wa— V fale oy u® 
—Rhiva ue Wi, , 
where Wa,=VE Vs. 


Then we take a covariant vector a, and a contravariant vector 07 such 
that ua, =v’a,=0, and multiply (6) by a v. On contracting, we get 


(7) | Reeve a 0" Op + Runes x VE 8; + Rue Qn VF 8, 
—Rrayw ay ue Wi—R “hae Op 04 Wi—Rtwa @,u® Wi=0 


where Wi= Wis. These equations have the form 


ROEDER SR M e 
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(8) Lye Op + Myo 83—Muy 85 
—CyoWi—Duw Wi+Du, Wo=0 
with Lyw=Mww—Moy,  Cyw=Dyo— Day , 
if we put De Rs Miuw=R pe dav® , 
Cre Rive Gath’, Dia sRisae Gat? , 


for we have Riyy.3=0 and Rioa=0. Then owing to the lemma proved in 
the appendix of this paper, we can conclude 
(9) MD 

These equations mean 
(10) RB daVP=AR upy Aa UF , 
in which @ is an unknown quantity and the vectors w, v2, and a, must 
satisfy | V2.e*|3=0 and u’a,=v’a;=0. But as the former inequality is 


satisfied by an arbitrary vector w* as long as it is sufficiently near some fixed 
vectors, only the latter equalities are to be taken into account. As we get 


from (10) Rgy Za VB’ =0, we have 
(11) Rrupy Via =u A, +B, , 


where A, and B, are unknown vectors, which will depend upon z#* and 0. 

Now we can take w* and v* arbitrarily, for we can take the equations 
Ua, —0, v’a,=0 as the equations confining a, but not w* and v*. If we put 
ur=S}, v'=83, we get from (11 Rl23=02, Similarly, we get Ri,;,=0 for 
i+=j, ixck, and this is true with almost any coordinate system? If we 
change this, we get, for example, R%,syABB}=0, where A=dy/ax, B? 
=0x/oy4. Put At=d)+16} 82, Bi=84—t8} 62, at the point (x**)) in question, 
then we get R%,y(544+ 284) (8; —26/ )83=0, that is, R2.23=R'13. Similarly, we 
get 


Oh ee oe = 
Roi =Regi= Oe =", 
RI = = te eee ee =) Ped 
petal sy oun? 


eeorsesee = eee eeesee 


and from these facts we get at last 
(12) 2 Rests Ae 8 A ° 
As this affinor must satisfy Ritwel=0, we have Ay=Ayyp. 


Substituting (12) into (10), these become Ays v? =a Aye A. In these equa- 
tions w* and v* are arbitrary, while A,y are independent of w* and 0’, for 


(2) It may happen that 2 =a make | VA. ur |=0, but remark that we can take such 
a coordinate system that | Vin u*|=-0 is fulfilled by any one of w= M5 ' =. 
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we have (12). Then we know that the affinor A,, has the form A,y=A,8B,, 
and as Ay, is symmetric, we get A,,=RB,B,. Hence R>.. has the form 
(13) ReERD; (8% Bu ー84 B,) . 


We postpone further discussions a little and treat 


§3. The case where | Vi * |=0 for any *. 


In this case we can find out a vector W, for each vector z#*, such that 
we have 


(14) cp UH" W=0. 

Then multiplying (5) by uw* W, and contracting, we get 

(15) Riise Wr O=Rrevo US Wy + Rne UX W, + Reine U We. 
Multiplying these by an arbitrary vector a“ and contracting, we get 
(16) faite ux Wa Cau W,+R ux Wa —R*, ux W, > 


where R>,=R>%,,. a". Multiply (15) by wk and contract, substitute (16) into 
the resulting equations. then we get 


CR U* Wr—R>. u* Wa u®) W, 
=(w Ro, u* W,— Ro, u* Ws 48) Wa. 
From these equations we obtain 
(17) wR, UX Wr—Rry u* Ws uP=aW,.. 


Multiply these by w’ and contract, then we get a*W, u’=0 because of Row) 
=0. This means aw*=0 or W,u’=0, hence we get 


Foy 

(18) u Wan u 

or 

(19) W, w=0, W, Rr, uz=aWy. 


If (18) is true, we have R%,,, a" u’=u7,, where 7, is an unknown 
vector. As ar is arbitrary, we get Rive w= Tew With unknown 7. If 
we put ~\=d}, we get Ri,;.—0 for 74¢1. We can obtain a set of similar 


equations, and as Riae)=0, we can conclude R>,,.=0. Hence this case must 
be discarded. 


Let us consider the case where (19) holds. 
Then multiplying (15) by u” and contracting, we get 


A 
Rep 1% uP=L- , 


because’ of W;, w=0 and Rou.” w=0. As wu are arbitrary, we get, as a 


ee on. ge ately 
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result of some calculation usual in tensor calculus, and such as used in 
obtaining (12), 


(20) Diet Rina Panter Baw dh ; 


where B,, is an unknown affinor. 


As we have Rwel=0, we get from (20) 
OR aise RO OD 
Exchange v and » and multiply the resulting equations by 2 and add them 
to the former ones, then we get 
SRrouv=(2Bro + Bay dp +2By 85+ Bon & 
But as Roe)=0, we get 


(Byot Boy) + (By + Buy 8+ (Buot+Bow&=0, 


from which we can conclude B,,+B,,—0. Hence we have 
1 
(21) Ree ney (Buy oa Bre 83+ 2By. 8) ° 


It seems that we can draw no more useful results from (5). Therefore 

we substitute (21) into the Bianchi’s identity Rd.ry.-0;=0, and get 

6 Sn Bry sot 8y (BB 0) 

+8 (Buo':v— Buy; 0) +8> (Buy; 0 Bue; 5) =O, 
Contracting with respect to X and o and taking into account of B,,+B,,=0, 
we get 
(22) 2 Byosut 2 (By; o—Byo;v=0, 
and taking the symmetric part with respect to pp and v, we get Buo:yt+ Bo: 
=0 for ~=>3. Hence By: is skew-symmetric in all indices and we obtain 
Buy; 0=0 

from (22). 

Differentiating covariantly, we get Buy:o;0=0, hence R%uc Bavt Rowe Bua 
=0. Substituting (21) into these equations, we get, after some calculation, 
2Buy Bort Boy Buc—Bov Buw=0. Taking the symmetric part of these equations 
with respect to v and o, we get Buy Bort Buc Buy=0, hence B,, must have the 
form B,y=A, By, with unknown A, and B,. But we know Buy=—By,, and 
consequently we obtain Bs=0 and Re=0. Then we get the 

Lemma: An affinely connected space A, with non-vanishing curvature 
affinor can not have the curvature affinor satisfying (21) for n=3. 

Thus we must discard the case | V2, w* |=0. . 


But here we have an interesting remark. 
Suppose that we have a group of affine motions with maximum. order 
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yr, nz<r<mn+n. Then we have equations (3) with 7=1,…, 2+%—?, but we 
can fill the vacant ones, i=n2+n—rt+1,-m, formally by 0xé&,+0xé&*=0. 
The determinant | Vi, w*| thus obtained is always null and we get the 

Theorem: An affinely connected space can not have a complete group of 
affine motions of order r, n2?<r<n?+n (Egorov, 1]. 


§4. Relations obtained from XR nic=0vend Xv: =O. 
We return to the equations (13) and rewrite them as 
(23) Rive=eAi (85 Ao —8. Ay), 
where e=+1. Substituting these into XR%,,.=0, we get 
[(XAy Ay + Au XAuI—[(XA, )Ay+ Ay XA] 83=0, 


hence 

(24) (XA, YAy+ An XA,=0 ’ 

from which we soon obtain XA,=gA,, and substituting these into (24) we 
find 


(25) XA, =0, 
that is, Ay;>o &£° +A &7;,=0 or 
(26) E87 OL Act+E* Au;.=0. 


As there are just ” linearly independent equations in (26), we know that 
(26) is equivalent to (3). 
Substituting (23) into X(Ri..:-)=0, we get 
XT Au;0( 8) Aw —83 Ay) + Ayu (8) Aw —83 Aye) =0, 
from which we obtain 


X(Ay;0) Aw ae A, X(Aa;0)=0 > 


hence 
X(Ay;)=0, 
which mean Ay:v;0 £%+Ao;y Ep + Ap;o £7 y=0 or 
(27) E57 CAk;y Sp t+ Aus ne JAE" Ausy;e=0. 


But these equations must be satisfied whenever our condition (3) or (26) is 
satisfied, hence we have 


(28) Pie: }.On T Aare Oy One Op aieeeicae. whey 
(29) Aas Lin flack ’ 


which mean that the equations in (27) are linear combinations of the equa- 
tions in (26). 


Contracting with respect to r and p, we get from (28) 
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(30) NAg y+ Ay e=Lay" Ag « 

Similarly, we get Ac:yt+7Au:.=Lii"A,, hence 

(31) Ag ;y+ NAy = Ag . 
From (30) and (31) we get 


(w= I) An y= ALay line YAg 
and 
(#2—1)Ay, RD — Lay”) Ak , 


and from these equations we can deduce 


(32) Ap: v=aAy Ay ° 
Substituting (32) into (28), we get 
(33) Lj) =a A, EA 


and substituting these equations and (32) into the right hand side of (29), 


we get 
Ax; Vv; c= 207A, Ay Ax . 


On the other hand, we have Ay;y:”«=2@?Ay Ay Ac+ @;x Ay Ay from (32). There- 
fore we get a..=0 and 


(34) a=constant. 
§5. The expression of an A, admitting a group of affine motions of order 
n2 in tensor form. 
All that we obtained hitherto is the set of relations 
Rww=eA, (83 Ay —83 Ay) 
(35) Au;v=@A, Ay, 
a=constant, Ay*<0, 


as the necessary condition for an A, to have a group of affine motions of 
order m2. The equations for the vector &* are 


(36) Apt aA, &* Ae =0. 


That (35) is the necessary and sufficient condition is found at once. Because 
of (35) we have 


pee ef =(m+1)! am cA, (83 Ay —83 Ay) Ao, Aon» 
and as we have XA, =0 already, we soon find } 
XR Saves i “300. 


The last thing to do is to prove the existence of an Ax satisfying (35), 
but this problem was already solved by Egorov by finding an example 
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[Egorov, 1]. I do not represent it here, but instead, in the next section, I 
solve (35) using a suitable coordinate system. As a result we have the. 

Theorem: An affinely connected space An with non-vanishing curvature 
affinor has no group of affine motions of order more than n2. It has a complete 
group of affine motions of order n? if and only if the condition (35) is satisfied, 
and such spaces really exist. 


§6. The property of an A» fulfilling (35) and the complete group of affine 
motions in this An. 
We find that A, is a gradient vector, A,=f,, for Ay,y—A,,,=0. Then 
we can take a coordinate system such that f=x”", and get 
(37) Ag=orn 
Because of (32) we soon find 


8) ; lO 


for i, j, k,--=1,-, n—1. 
From the expression of Ri we get, for \=i, 


(39) Penn one —U Gn ek Ge an=0, 
(40) D5 deh ih = 
(4) nA nel ab EL 
(42) Eno! wk Lael ok mal wht gil oneoks 


Take a hypersurface x%=c and name it H(c). In each Hc) we can find 
out a set of m—1 independent parallel vector fields because of (39), and this 
means that we can find out a coordinate system such that 


throughout the Ay. Besides, as we can choose the origin, #=0, in each H(c) 
such that the locus of this point for variable c forms a geodesic in the An, 
we can put 

(44) ri,=0 for  x=0, 


because of + Ii, x" P= 

In taking a coordinate system, we have another freedom, the direction 
of the (n—1)-ple in, H(c), formed by the vectors uj=Si,---, u/,.;=84-, at the 
origin xi=0. We can put them parallel along the geodesic xi=0. Then we 
get 


(45) Ti,=0 for -xe=0, 
because of 4,,,+/74,4%=0 and u)=d}. 
Substituting (43) into (40) and taking into account of (45), we get 


ee a 


a ee 
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(46) P05. r= 0 
throughout the Ay. Similarly we get from (42) and (44) 
(47) Ti ,=—exi. 


We get nothing more from (41). ’ 

The last parameter 1%, is calculated from (32) when we put p=v=n, 
2 — a. 

All that we found is that the parameters J"), are given by 
(48) I't,,= —exi, Ly —-a other 73,0 


We get no more relation from (35), and (48) is the complete solution of (35). 


Now we are ready to find out the complete group of affine motions in 
the A» fulfilling (35). 


From (36) we get £%,+a£"6=0, which becomes, because of (48), €” ;=0, 
&" ,=0, hence 
(49) &"=constant , 


while from (1) we get 


(50) SEN EL tye! enw pe Cy 0. 
If we put A=z in (50), we get 

E* 5, ,=0 

と k=O 


Ei nt Dina b°—Ln Ea =0, 
hence 
Bj Kj wl fF HIGH EP OD, 
fi’ + af’ —efi=0. 
If we put A=~z in (50), we get no more relation. Hence we get 
61) | &i= Ki 4i+Aiexp (a, x”) + Bi exp (a x") 
g"=K, 


or 


| £i= Ki xi+(Ai+ Bi x) exp (atx") 
(52) 


n=K, 


where Ki, Ai, Bi, and K are 2 arbitrary constants, and a’s satisfy 1°+ aa, 
—e=0, «2+ aat2—e=0, Ai A2; a? + aAa—e=0, 2a+ a=0. 
A finite transformation of the group is easily found to be 


xi=Ai xi + Ai exp (a, x”) + Bi exp (a2 x") 
(53) 


gli Aa “ - 


10 Y. Muto 
or 
xix Al. xi + (Ai+ Bi x”) exp (Car) 
(54) | ) 
rn— M+ A, 
with | 4 | >0. 
We can resume the results in the 
Theorem: Az affinely connected space An with non-vanishing curvature 


affinor admits a complete group of affine motions of maximum order if and 
only if the equations 


Rouveo = eAy (8 Aw — $s A,) > A,; y=aA, A, 5 a=constant 


are satisfied. Then the order is n?, and we can find out a coordinate system 
such that the coefficients of connection satisfy 


Pi ee oy Other fg 0. 


and a finite expression of the group is given by (53) or (54) according as the 
roots a and as of the equation x?+ax—e=0 satisfy a,*= a2 or A|=A2=a. 


APPENDIX 


We prove the 
Lemma: Jf we have the tensor equations 


(Al) Lyd + Myo 8 — Muy 83 


—Cyo We — Duo WS + Duy W2=0 
and 
(A2) Lyo=Myo—Mun, Cyo=Dyo—Day 


for n= 3, then we can conclude 


(A3) Lins D0 
Or 
(A4) Muy=aAD,y ’ 


where a ts a scalar. 
Let us begin with the case m= 4, for in this case we can prove the 
lemma easier. 


First we assume D,,=0. Then we have C,,=0 and (Al) becomes LZ, 5% 
+ My 8) — My, &$=0, from which we can find that (A3) is true. 
We next assume D,,+D,,=0. Then we have 


© Ly 88+ Mi 88— Myy 88 —2D,4 W2— Dy, Wi + Duy W2=0. 
If we have some vector Wr satisfying M,, Wr W’ 3=0, then we get 
Muy Wr WY 8= We (Lyn + Myo) W’—3 CW WD,, W), 


Eg A EE gy 


a 
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hence 3=WT, + WV, with some vectors 7, U*, and V,,, which is not true 
for n=3. Therefore we know that Mj+M,,=0, and (A1) becomes 


2M,, $+ Mus 8) — My, 83—2D,. W2— Dy Wi+ D,, W=0. 
Taking the symmetric part with respect to w and v, we get 
Mi 8 十 M,,.. d= Dy» W} — Dro Ws =0 ’ 


and multiplying these equations by some vectors c’ and d® and contracting, 
we find that the affinor W* has the form 
(A5) Wi=adi+ ST, + UV,. 

If we assume D,y+ Dy.3<0, we can find out some vector WA such as. 
Diy We WY *=0. Multiplying (Al) by W“W” and contracting, we get (A5) 
again with some S‘, 7,, U*, and V,. 

Substituting (A5) into (Al) and putting 1,,=L,.—aCyo, Myvy= Myuy—AD ays. 
we get 


(A6) Dw Oi Myce Oy — My 5 十 […]S*+[…] N= 
If we multiply (A6) by a vector a, satisfying $a,=U*a,=0, we get 
(A7) ba dp st Mico a;— Muy aw =) . 


Take such a coordinate system that at the point (%*)) in question we have 
a=6!. Then if we put o=1, Bl, v7=1, we get mjr=0 for j+-1, k3¢1, 
while if we put w=1, w=1, v-—1, we get 1;,;—m,;=0, hence mj,=2m,; for 
jx<l. As n=4 we can find out another vector 6, satisfying $b,= Wb,=0, 
and let the coordinate system be such that a=6h, b,=8, at (oD, Then 
we can deduce m,,=0, which mean (A4). 

If n=3 we can not find the vector 0, mentioned above. But according 
to a well-known property of matrices we can choose such a coordinate system 
that the components of the affinor Ws satisfy one of the following sets of 
equations ; 


1) Wi=ao., .Wi=a's. We=a''d., 
aaa ae els, 
CII) Wi=ad, W=ae, 
(III) a= a= a Wig + als, 
: asa’, 
(IV) Wi=adi, Wi=dita&%, Wp=8%+a8, 
at (x). 


Let A, uw, V, and @ take on various values out of 1, 2,3. We obtain some 
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sets of linear equations with M,, and D,, as unknowns, and it is a matter 
of simple exercises to solve them. We find 


(1), CID M0 ’ 
dD, (Iv) Miw=aDyr ’ 
hence (A3) or (A4) is true. 


In concluding the paper, the author wishes to express his hearty thanks 
to Professor K. Yano, whose kind encouragement and valuable suggestions 
have made the completion of the present paper possible. 
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I could not read the first two papers, which form a part of Egorov’s excellent disser- 
tations. According to the Math. Rev., Vol. 14, he obtained among other conditions (b) 
Rap=(1—7) 028, (C) @:8)=Cw28, Where 2.«;p Satisfy 2«#;ra2y7=0. But then we have 
a; By =eg: yAB, hence #;s=yarg with some vector vs. (c) becomes vers typ ra = 
2C)e28, from which we get we =Cle¢ , hence 2¢ ; 8=Cla 18. Perhaps these equations cor- 
respond to my (32). 


Added im proof: According to the Math. Rev., Vol. 14, No. 11, G. Vranceanu studied 
the same problem, but I could not read the original paper. 
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Preliminary Report on the Magneto-thermoelectric 
Power of Nickel 


By 


Nahonori MIYATA and Zenya FUNATOGAWA 


Synopsis 


The magneto-thermoelectric powers of nickel are measured by the same 
method as in the case of iron. The results show positive signs and that (1) 
the magneto-thermoelectric power has a maximum at about -90°C and 
then decreases with decreasing temperatures and may tend to zero at 
0°K, and (2) the magneto-thermoelectric powers between -90°C and about 
300°C decrease with increasing temperatures but (3) in the neighbourhood of 
the Curie point they are negative signs and are produced by the additive 
magnetizations, which are proportional to the external magnetic fields, ex- 
pected from the Weiss’ theory. 


§1. Introduction 


There are many investigations concerning the changes of thermo-electric 
forces of ferromagnetic substances such as iron and nickel, produced by 
magnetizations since the phenomena were discovered by W. Thomson. Re- 
cently one of the authors investigated the changes of thermoelectric powers 
of iron single crystals due to magnetizations and obtained many interesting 
results. We are now investigating the phenomena in the case of nickel by 
using the same apparatus and method described in the previous papers.‘? 

The magneto-thermoelectric force is the thermal e.m.f. of the thermocouple 
consisting of a demagnetized specimen and the magnetized one, and isa 
difference between the two electromotive forces, each of which is measured 
against copper respectively, first for the demagnetized state of a specimen 
and then for the magnetized state of the same specimen subjected to external 
magnetic fields under the same temperature gradient in the specimen. In 
many previous papers the effects have been measured under the condition 
of the wide temperature differences between the cold and hot junctions, 1. e., 
100°C or more because the electromotive force of this couple is very small. 
It is unreasonable to discuss such results from theoretical view points because 
the magnetizations and thermoelectric powers may not change always with 


temperatures in a linear relation, so that we must define and measure the 


apt AT AT 
thermal e.m.f. as mean values between the temperatures + mer and T— i 
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as in the case of ordinary thermoelectric powers. For simplicity, such defined 
e.m.f. is called the “ magneto-thermoelectric power” or the “ magneto-thermo- 
electric force per degree” at the given temperature 7. 

It is the purpose of this experiment to know the temperature dependency 
of the magneto-thermoelectric power of nickel and to compare it with the 
results of iron and also to give some fundamental knowledge to clear the 
mechanism of the interaction between the 3d- and 4s-electrons. 


§2. Samples and results of measurements 


The specimens of nickel are cylindrical rods, their diameters 1.00mm. 
and lengths 25.3 cm. 


The magneto-thermoelectric power is defined by ee a at T, as men- 


mag 


tioned above, where 4£ noted 
by a galvanometer deflection 
is the induced thermoelectric 
force by the magnetization of 
the specimen when the tempera- 
tures of each end of the speci- 


men keeps T+ Ar and T— 


respectively, under the condi- 
tion that the ordinary thermal 


e.m.f. of the couple consisting 
of the specimen and lead wires 
is conpensated by a potentio- 
metric connection. We have 
measured the temperature 
range between the boiling point 
of liquid oxygen and one above 
the Curie point under the con- 
dition that the 47 is about 
5°C~35°C and the maximum 
longitudinal external magnetic 
field is about 700 Oerstead. 

A typical example of the 
results is shown in Fig. 1: the 
Fig. la is the magneto-thermo- 
electric power, which is positive 


Fig. 1. (a) shows the magneto-thermoelec- 


: tric powers with the magnetization and (b) 
magnetized metal to de- shows the magnetization curves. 


sign, i.e., current flows from 
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Fig. 2. The absolute thermoelectric power of nickel. 
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Fig. 3. Experimental relations between the magneto-thermoelectric 
powers and the magnetic fields at various temperatures. 


‘magnetized one at cold junctions and the Fig. 1b is the magnetization curves 
of the specimen. We can not distinguish the differences among the kinds of 
the demagnetizing methods such as A.C. demagnetizing- and thermal ‘de- 
‘magnetizing- and other methods. The absolute thermoelectric powers of 
nickel on the temperature range between the boiling point of liquid oxygen 
and about 550°C are shown in Fig. 2, results of which are always negative. ' 
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Fig. 3 shows the functional relations between the magneto-thermoelectric 
powers on various temperatures and external magnetic fields. Fig. 4 shows 
the nearly saturated values quoted from the Fig. 3, where some corrections, 
as will be described immediately, are given to the values near the Curie 
point. We know from these results that (1) the magneto-thermoelectric 
power has a maximum value at about ~90°C and decreases with decreasing 
temperatures and it may be true to tend to zero at the absolute temperature 
zero, and (2) there is a linear decreasing relation with increasing tempera- 
tures between the maximum point (-90°C) and near the Curie point, but (3) 
there appears a sudden decrease and shows a minimum negative value at 
the Curie point. 


§3. Discussion 


As pointed out in the previous papers the saturation values of magneto- 
thermoelectric powers are indefinite and are affected by the distributions of 
magnetic domains in a demagnetized state of the specimen. In our case, 
however, we may assume that there is no difference because the distributions 
of magnetic domains in demagnetized states are not so much affected, as 
described above, by the kind of demagnetizing methods on account of the 
some residual internal stress and comparatively small anisotropic energy 
constants. 

We observe that the maximum saturation value of magneto-thermoelectric 
power appears at about -90°C as one of the authors has observed in the 
case of iron. It is reasonable that the maximum may appear at a certain 
temperature because the magneto-thermoelectric power must tend to zero at 
0°K as was proved in the previous papers from the stand-point of the 
Nernst’s heat theorem which states that the thermal e.m.f. must be zero 
at 0°K. 

In the neighbourhood of the Curie point, since the additive magnetization 
proportional to the external field becomes more conspicuous than the magneti- 
zation due to the change of distribution of magnetic domains, the related 
phenomena such as magneto-resistance effect and magneto-caloric effect etc. 
must be distinguished from the effects which appear far below the Curie 
point. In our case, because the temperature gradient in the specimen may 
be changed by the magneto-caloric effect, the observed values of the magneto- 
thermoelectric power must be corrected by the measurement of temperature 
variation produced by this effect. Thus, for example, the linear part of 
the increasing of the magneto-thermoelectric power with the increasing 
magnetic fields as shown on the curve T=380.9°C 4T7=22.0°C in the Fig. 3 
can be explained as due entirely to the magneto-caloric effect from our 
measurement and some calculations. The relation between such corrected 
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values at H=570 Oerstead of the magneto-thermoelectric power and tem- 
peratures is shown in Fig. 4. Nevertheless, a negative deep valley appears 


4E 
Cp) mag. 
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Fig. 4. The saturation values of the magneto-thermoelectric powers vis. temperature. 
The saturation values at low temperatures are obtained from the curves in the Fig. 3 
by the extrapolation method and those in the neighbourhood of the Curie point are 
the corrected values at H=570 Oerstead. 


in the neighbourhood of the Curie point, produced by the increasing additive 
magnetization expected from the Weiss’ theory. 
mus Fig. 5 shows in more detail the 
(C.R) ; magneto-thermoelectric powers in the 
0 A OO 100 EE neighbourhood of the Curie point with 
the magneto-caloric effect measured by 
us (0o) and also with the magneto-resist- 
ance effect measured by Gerlach? (x) 
in suitable scales in order to coincide 
with the magneto-thermoelectric power 
at the Curie point. We know that from 
the coincidence of the effects these are 
caused by the similar origin, probably as 
Fig. 5. The magneto-thermoelec- mentioned above by the additive magneti- 
tric powers in the neighbourhood of zations proportional to the magnetic 
the Curie point (C.P.). x and © fields, and also that the decreasing of the 
anpthe jtigure show \the | magneto- absolute thermo-electric power by the 


resistance effect and magneto.caloric Son wee 
effect, respectively, indicated in sui- appearance of magnetic domains below 


table scales in order to coincide with the Curie point as shown in the Fig. 2 
the magneto-thermoelectric power at corresponds to the decreasing of the 
the Curie point. magneto-thermoelectric power produced 


by the additive magnetization in the neighbourhood of the Curie point. If 
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we put that SR and 8E are the change of the electric resistance and the 
thermoelectric power, respectively, by the appearance of spontaneous magneti- 
zation at 20°C just below the Curie point and 4R and 4E are the change of 
the corresponding quantities by H=400 Oerstead at the Curie point, then the 
next experimental values are obtained 


AR 4E 
— -2 —— -2 
SR =13X10"7, SE =1.2x10-2. 


§4, Summary. 


We summarize our experimental results as follows: 

(1) The magneto-thermoelectric powers far below the Curie point are 
produced by the magnetization due to the change of distribution of magnetic 
domains and those in the neighbourhood of the Curie point are produced by 
the additive increasing magnetization proportional to the external field. 

(2) At low temperatures both the magnetization and the magneto- 
thermoelectric power can not be easily saturated on account of the increasing 
of magnetic anisotropic energy constants. 

(3) The magneto-thermoelectric power has the maximum value at about 
-90°C and then may tend to zero at 0°K. The linear relation between 0°C 
and 300°C coincides with that measured by Broili®. 

We are now measuring the magneto-thermoelectric power and magneto- 
striction of single crystal of nickel. The results, which are very interesting, 
and more detail discussions will be published recently in the Jour. Phys. Soc. 
of Japan. 
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Temperature Dependency of Magnetic 
Anisotropy of Cd;Mg 


By 
Nahonori MIYATA 


Alloys of Cd and Mg show intimate relations between their order-disorder 
transformations and magnetic properties as reported in the previous paper. 
It is interesting for us to investigate their anisotropic characteristics by 
using the single crystals, since they have the hexagonal structures. 

The experimental method was similar to what had been used by Rao 
and Sriraman for Cd and Zn.” A few results are shown in the figures, 
where X, and X- are the atomic susceptibilities along a- and c— axis of the 
crystals in the equilibrium states respectively and %’s mean the atomic 
concentrations of Mg. The anisotropic changes occured with the construction 
of the superlattices. For alloys which had excess Mg than Cd3Mg the two- 
steps changes appeared along both principal axes, but their features were 
quite different. 

The temperature dependency of the mean values Xsoy= 4 (2X%a+Xc), which 
should give the character of the polycrystals, was in agreement with the 
previous investigation.) 
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Study on Lichtenberg’s Figures by Means 
of Color Films. (3rd Report). 


By 


Bunji ARAI 


Synopsis 
The background of Lichtenberg’s figures obtained on color films was not 
colored black, but showed some color. What acted on the color film? The 
problem was investigated by the present writer, and it ‘was found that these 
phenomena originated from the characteristics of color film of this country. 

- Some red streamers appeared often at tips of negative figures by color 
film. The writer considered that this fact would be due to the perpendicular 
construction on Lichtenberg’s figures. To confirm this fact, experimentally, 
the writer obtained micro-photographs of section of streamer in Lichtenberg’s 
figures by ordinary black and white film, and observed the distribution of 
silver particles in, an emulsion layer at any parts. 


1. Introduction 


The background of Lichtenberg’s figures obtained on ordinary films or 
dry plates by orthochromatic emulsion had no color and was transparent. 
But the background of Lichtenberg’s flgures obtained on color films by 
applying impulse, turned red over the whole film surface, as shown in 
Photo. 1. About this fact, the writer described already in former reports.) 
Now, even if the red fog is shut out by the black paper by using tungsten 
type, yet the background of figures appeared in dark olive green, as shown 
in Photo. 2-a, and by using sun light type dark violet, as shown in Photo. 
2-b, which naturally should be black. The writer considered that the phe- 
nomenon was attributed to either of the action of strong electric field or 
character of color film. The judgement was made by the method to be 
described in next section. 

The tips of negative figures obtained on color films in an ordinary posi- 
tion turned out usually red streamers. This fact seems to show that a red 
sensitive emulsion layer in the nearest celluloid base was excited by charged 
particles, as described in former reports.):2) Then, he supposed that also in 
_the case of. black and white, film, the silver particles at tips of a negative 
figure are in the nearest layer to celluloid base, as shown in Fig. 1... 

Previously, C. E. Magnusson?) observed a slender prong of Lichtenberg’s 
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figures Ly the microscope in perpendicu- 
lar direction of surface of the dry plates. 
He described that “ Visual examination 
of both the positive and negative figures 
by means of a compound microscope is 
more satisfactory, and that to study the 
microphotographs as for the higher 
magnification the thichness of the emul- 


Fig. 1. Supposed diagram of sec- sion film is greater than can be observed 
tion at tip of negative Lichtenberg’s without the change of focus. The dark 
figure. h specks are distributed through the emul- 
No. Name Thickness = fil d 5 

Pi Alin Du sion film and not merely in a surface 
2 Sensitive emulsion layer 20 p. layer”. The author researched directly 
3 Celluloid base 100 p. the distribution of silver particles from 
4 Antihalation layer 15 yp. the cross section of the sensitized film 


differing from Magnusson. 


2. Experimental Apparatus for the Discrimination 
of Color in the Background 


The connection for supplying the impulse to the films set in these cameras 
is shown in Fig. 2 | 

Through the experiments, the sphere spark-gap G; (25cm ®? brass ball) 
was used at 0.95 cm (30kV peak value, 25°C, 760mm Hg.), the impulse 
applied by G; was chopped by the sphere gap Gz». Gs; (2.0cm¢ brass ball) 
was used at 0.69cm (21.0kV peak value, 20°C, 760 mm Hg.). 

These experiments were tried by the film of tungsten type only, because 
the cut film of sum light type was not prepared in market. For these ex- 
periments, Fuji Color Cui Film (tungsten type, non halation, safety base, 
Em. No. ATC-190-G) was used. The development of these color films was 
submitted to the Research Laboratory of the Fuji Photo-Film Co. Ltd. 

For the glass plate in the special camera, the glass plate of dry plate 
was used. The glass plate was a cleared one which removed the emulsion. 
The thickness of this plate was 1.20 mm. 

The metal plate in the special camera which was described in Fig. 2 of 
the former reports? was cut out from the half part, as shown in Fig. 3. 
There, the distribution of strength of electric field was varied from left to 
right. The colors of the parts left and right in the background of the ‘HEUER 

obtained by the cameras were compared. _ 

Also, the camera with the electrodes of point-ring arrangement was 
' prepared, as shown in Fig. 4. The color inside the ring of the ‘figures 
obtained in this camera was compared with that outside the ring. = 
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——§fo D.C. high - tension source. 
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Fig. 2. Schematic connection diagram of impulse circuit to the special camera. 


Sphere spark-gap (2.5 cm ¢ brass ball, gap was fixed at 0.95cm, 25°C, 760mm Hg. 
30kV max.). 

Sphere spark-gap for chopper of impulsive wave (2.0cm ¢ brass ball), G2< Gj. 
Large Leyden Jar (0.0038 pF x 2). 

Liquid resistance for damper of impulse (solution of CuSO4: 0.4KQ). 

Fixed resistance (HOKUOHM-—resistance ready-made as a radio parts, 1MQ~ 10, 
+ watt type). 

Black paper (cover of color film, 0.09 mm thick). 

A group of sensitive emulsion layers of the color film. 

Glass plate (cleared by removing the emulsion of aged dry plate. 1.20 mm thick). 


“ ‘To D.C. high tension source” leads to Gaiffe-Gallot et Pillon’s apparatus.) 


“Fig. 3. Details of special camera Pdesumed for this experiment. Metal plate 
(P.E.) was cut out from the half part. , 


‘This camera was made by using an empty box of Agfa’s Isochrom dry plates. : 
The. whole stface of this camera was coated with bakelite varnish for keeping perfect 


insulation. 
I : Bakelite plata (0.26 cm thick), B-: Black paper, - mo td 
L : Bakelite lock nut. (1.0 cm $), _ Ee : Glass plate, yn 


"PQ: Copper threaded rod (0. 20 cm #), J.X.Y: Paper tox. (size—10.8 cm x 


Fi: Color Film, _. 83cemx19cm), 


C.E. Brass cylindrical electrode © 70 Th: Terminal "for high tension side. ' 
cm $) Threaded in P. 2 i 9: Te:Terminal.for earth side: Doni ns 
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Fig. 4. Details of special camera 
with the electrodes of point-ring 
arrangement. The notation is the 
same at that for Fig. 3. 

But, I: Ebonite plate (8.3 cm x 
10.8 cm x 0.5 cm), 

PQ: Brass rod electrode (0.30 
cm ¢), Threaded in I and locked 


3. Results and Discussion for the Color 
of the Background 


Photo. 3-a and 3-b show respectively 
the positive and negative figures obtained 
as shown in Fig. 3. In this case, no dif- 
ference was found between the colors of 
left and right parts. Photo. 4-a and 4-b ° 
show the figures obtained by the cameras 
shown in Fig. 4, respectively. Photo. 4-a 
show the figure in which the center point 
was connected to the plus and the outside 
ring was connected to the minus (earth 
side). Photo. 4-b gives the inverse con- 
nection. In this case also, no difference 


by L. 
Ri: Brass ring electrode (Inside 
dia.—5.0 cm ¢, width—0.50 cm). 


in color was observed between the inside 
and outside. 

In order to make sure this phenomenon, 
the writer observed the color of the film to which no impulse was supplied. 
There, no difference was recognized between the applied film and the not- 
applied one in the film of the same package. They were developed at the 
same time. 

Then, the writer recognized that the colors of the background of Lichten- 
berg’s figures which appear on the color films are the character of the color 
film in this country. 

As seen in Fig. 7-a and 7-b, orange color appeared at the boundary of 
the metal plate. About this phenomenon, the writer thought that the red 
sensitive emulsion layer (deepest and nearest layer to celluloid base) and a 
part of green sensitive emulsion layer were excited by the strong electric 
field. 


4 _ Sample and Praparat Making of the Section of» 
Lichtenberg’s Figures in Film 


At first, the cut film (Fuji Portrait Pan.) was tried to be cut. But, it 
was difficult to cut a thin piece at uniform thickness. Its celluloid base was 
too thick (about 200) compared with an emulsion layer (20). Then, for 
this research, Roll Film (Fuji Neochrom Brownie, orthochromatic, F 120, non 
halation, Em. No. 021-D) was used, its celluloid base having a half thickness 
of cut film. The orthochromatic film was taken by reason of no necessity 
of the’ cover of black paper! to protect the fog over the whole ‘film surface, 
and of easy dark room operation. This roll film was cut into about 10cm 
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length. The emulsion of this film was put in an upper position, and was 
stuck on the glass plate by photo-paste at either end. The glass plate was 
a cleared one which removed the emulsion from the dry plate. These glass 
plates were 1.20-1.23 mm in thickness. After natural drying of the photo- 
paste and until perfect sticking in the dark room, it was put in the camera 
(as shown in Fig. 2 in former reports2). In order to apply the same crest 
voltage, to obtain clear and sharp positive and negative ordinary figures, 
. and to apply the pulse at the same time on the films set in these cameras, 
the connection was arranged as shown in Fig. 5. 


fos te high- tension source. 
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7 
Fig. 5. Schematic diagram of the impulse circuit. Connection to these cameras, 
to apply the same crest impulse voltage, to obtain positive and negative figures 
at the same time. The notation is the same at that for Fig. 2. 


See Photo. 5-8. From the emulsion layer of the film which was put. in 
down position and stuck on the glass plate, unsharp inverese figures to 
positive and negative were obtained respectively as shown in Photo. 9, 10. 

The writer picked up a typical streamer in these figures, and cut out 
strips including a streamer along other streamers from these films. Some 
samples of thin pieces were prepared by the following two method. 

i) One of the method of ordinary botanical praparat making was taken. 
These strips were put in paraffin inbed (melting point: 52°C), and then the 
samples of thin pieces in the inbed were cut out as thinly as possible at the 
desired parts perpendicularly to the streamer by a safety razor blade with 
free hand. These thin pieces were put on a slide glass, and paraffin was 


removed by xylol. A temporary praparat was made by being soaked in. rh 
under a cover glass. 
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ii) The other was a extemporaneously devised by the writer. The strip 
was inserted between two thick celluloid plates (about 3mm thick), and its 
part desired to be cut was put at the end of the celluloid plate. These 
plates were pressed strongly and the part which bulged out of the strip was 
cut out by a safety razor blade with free hand. The end surface of the 
thick celluloid plate was used as the guide. Next, the strip was bulged out 
as the thickness of thin piece, and was cut in the same manner. These thin 
pieces were put on the slide glass and were soaked in a drop of water. The 
temporary praparat was preserved in this state, under the cover glass. 

The samples were taken at the vicinity of the electrode, the intermediate 
of streamer, the tip, and the outside of tip. 


5. Preliminary Discussion for the Sectional Distribution of the 
Specks in the Film of Lichtenberg’s Figures 


The numerical figures given at the side of microphotographs show the 
distance from the cylindrical electrode in mm. 

Photo. 5, 6 show the series of micro-photographs by method i). The 
sample of ordinary positive figure near the electrode has a very dense 
distribution of the specks as shown in Photo. 5-0, and in other samples these 
specks distributed in the upper layer of emulsion. In the ordinary negative 
figure, the same distribution as that of positive figure was observed near the 
electrode, as shown in Photo. 6-0. But in the tip, these specks distributed 
in the lower layer of emulsion nearest to the celluloid base. 

The series of micro-photographs by method ii) are shown in Photo. 7-10. 
Photo. 7, 8 show ordinary positive and negative figures, Photo. 9, 10 show 
inverse positive and negative figures respectively. 

The distribution of specks in the lower layer .near the celluloid base was 
observed clearly in the sample of the unsharp streamer in Photo. 7-11.5. 
This fact corresponds to the unsharp red streamer in the color film in former 
reports).2, But in Photo. 8, the tendency as seen in Photo. 6 was not ob- 
served. Because all tips have not been observed, the above description cannot 
be insisted. 

In any case, the color distribution of streamers obtained by the color 
films as shown in Fig. 4, 6, 8,10 of former reports?) are explained by the 
vertical distribution of the specks obtained here. In Photo. 9, 10 the samples 
near the electrode have the dense speck distribution over the whole ‘range 
from the lower to upper emulsion layer. In the tips, the specks distributed 
at the lower layer only far from the celluloid base. These coincide with the 
color figure as shown in Fig. 5, 7, 9, 11 of former reports” For comparison, 
Photo. 11 shows the micro-photographs of the section of films by ordinary 
light. 
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From the above observation, the writer considered that Lichtenberg’s 
figures have vertical construction. 

This report is a summary of the lecture made at the Special Meeting on 
Discharge Physics, held by the Physical Society of Japan, Oct., 16. 1953. 

In conclusion, the writer wishes to express his hearty thanks to Dr. 
U. NAKAYA, Dr. C. MAGONO, Dr. T. ASAHINA (National Science Museum) 
and Dr. S. FUJISAWA (Director of Fuji Photo-Film Research Laboratory) 
who made suggestions to him throughout this work. 

The writer is also indebted to Assistant Prof. S. NAGAMI (The Botanical 
Dep.) for praparat making, and to Mr. T. ONO (Tokyo Sibaura Electric 
Works Co. Ltd.) for his supply of materials etc. in this research. 


Oct., 29. 1953. 
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Photo. 10-+15. 
Micro-photographs show cross sections perpendiculaly to streamers. Numerical figures 
written at the side of microphoto. show the distance (in mm) from the cylindrical 
electrode. Scale of min. devision show 10 wp. 


Photo. 10, 11: 
Sample: Paraffin inbed, soaked in xylol, 
Microscope and Camera: YASHIMA (x 150)—Hand made camera (x75). 


Photo. 12-15: 
Sample: Soaked in water, 
Microscope and Camera: WINKEL (x 154)—Canon ( x 45). 


Data of lichtenberg’s figures shown in Fig. 10->15 
Data and Time: Sept. 23. 1953, 12h45m->14h45m. 
weather: Rain, Room Temp.: 22.6°C, 
Atom. Press, (mm Hg): 762.6>761.9, Humidity (%): 75-76. 
G1: 0.95 cm (30kV), Gz: 0.50 cm (18.0kV), Cs: 0.0038yF x2, r:04kQ, R: 107Q. 
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Photo. 5. Ordinary positive figure Cx1.5) and its microphotographs (x 175) 
at cross section of streamer (x 2.5). 
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Photo. 6. Ordinary negative figure (x 2.5) and its microphotographs ( x 175) 
at cross section of streamer. 
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Photo. 7. Ordinary positive figure and its microphotographs 
at cross section of streamer. 
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Photo. 8. Ordinary negative figure and its microphotographs 
at cross section of streamer. 
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"Photo. 9. -Inverse positive figure and its microphotographs 
at cross section of streamer, 
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Photo. 10. Inverse negative figure and its microphotographs 
at cross section of streamer. 


4 


: 


ull 


tie BS Fv ee re 


Photo. 11. Micro-photograph at section of sensitized film by ordinary light. 
Neumerical figures show the notation, they are the same at that for Fig. 1 
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Photo. 1.: 


Photo. 2-a.: 


Photo. 2-b.: 


Photo. 3-a, 


Photo. 4-a: 


Photo. 4-b: 


Positive and negative Lichtenberg’s figures (15.6kV) obtained on 
color films, whole film surface turned red. 

Lichtenberg’s figures (21.0kV) obtained on tungsten type color 
film, the background appears dark olive green. 

Sunlight type (14.8kV). Both figures shows inverse negative 
one. 
3-b: Positive (dendritic) and negative (tertially) figures (21.0 kV ) 
obtained as shown in Fig. 3. Compare the colors of the parts left 
and right in background of the figures. 

Center point was connected to the plus high-tension terminal 
(21.0 kV ), and outside ring was connected to the earth. 

Inverse connection of Photo. 4-a. Both figures obtained by the 
camera shown in Fig. 4. 
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Photo. 2-a. (—). 


Photo. 1. (+&-—), naked. 


Photo. 2-b. (=). 


Photo. 3-a. (+). Photo. 3-b. (—). 


Photo. 4-a. (P:+, R: —), naked. Photo. 4-b. (P: —, R: +), naked. 


Photo. by B. Arai 
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On the Falling Velocity of Solid 
Precipitation Elements 


By 


Choji MAGONO 


Abstract 


Falling velocity of various solid precipitation elements was discussed 
theoretically, and simple relations between the dimension, the density, and. 
the falling velocity were derived, utilizing well known aerodynamic principles. 
and constants. 


Introduction 


It is considered that almost all rain drops originate from solid precipita- 
tion elements, and the rates of growth of these elements depend on their 
falling velocities. So the falling velocity seems to be one of the most im- 
portant problems in the investigation of the mechanism of the precipitation 
phenomena. Since prof. Nakaya? observed the falling velocity of snow 
crystals, no discussion has been done for the falling velocity of solid precipi- 
tation. To calculate the rate of growth of snow precipitation, it is desirable 
to formulate the falling velocity of various precipitation elements by their 
dimensions. To discuss the-falling velocity theoretically, it is mecessary to 
know their. falling states and their densities during the fall. The author 
observed the solid precipitation elements in their falling states by a strobo-- 
scopic camera?,, but their densities were then not measured. During recent 
two seasons, many measurements of the density of various snow flakes were 
made by the photographic method?), and the theoretical discussion became 
possible. It is the purpose of this work to search the relation between the 
falling velocity of solid precipitation, utilizing the aerodynamic principles. 


§I. Falling velocity of snow crystal of plane type 


a) Dendritic snow crystal. 

In the previous report?, it is ascertained that the snow crystal of plane 
type falls keeping its plane horizontally. When the snow has reached its 
terminal velocity u, its weight is balanced by the air drag, 


sh(o—p)g=Cps $ pu’, (1) 
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where s is the area of its horizontal cross-section, # its thickness, o its 
density, p the density of air, g the gravitational acceleration, and Cy, is the 
drag coefficient of the circular plate. Thus 


- sy Sa 2ho—P) , (2) 
Cp p 

From Nakaya’s observation, the thickness of this type was nearly 0.0011 cm, 

being independent of the dimension. o is considered to be 0.9g cm-3. Cy is 

assumed here to be 1.2, (Reynolds’ number* 102~103). This assumption is 

consistent with Dr. Spilhaus’ paper... On the assumption, the terminal velocity 

of dendritic crystal of plane type was calculated from (2), 


u = 35cm sec"l, 


which indicates that the velocity of plane dendritic crystal is independent of 
its dimension as Nakaya had pointed out!. And this value agrees fairly 
with the observed data as shown in Fig. 1. The difference between the value 
observed and calculated may be due to the estimation of Cy», partially to the 
fluttering motion of the snow crystal in falling state. 
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* Strictly speaking, it is impossible to select. Cp unless the falling velocity is known, 
‘because the drag coefficient varies with Reynolds’ number. In this paper, the falling 


velocity of snow crystals requisite for calculating Reynolds’ number were temporarily 
assumed to be 1m sec-~l, 
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b) Rimed crystal. 


Because no observation has been made for the density of rimed crystal,. 
it is difficult to estimate the thickness and the density requisite for calculating 
formula (2). Nakaya gave following formulae for individual mass of crystal, 


for plane dendritic crystal m = 0.00038 a3, (3) 
for rimed crystal m = 0.0027 d2, (3!) 


where m is the mass of the crystals in g, d the diameter of the crystal in 
cm. The author considered that the thickness of the rimed crystal observed 
by Nakaya was proportional to its mass, therefore, using (3) and (3 & is 
given, 

0.00038 d2 


= 0.0079 cm. 


In this case also, o is equalized to 0.9 g cm-3, for the thickness calculated 
here is equivalent to a compact ice plate. Thus the falling velocity of rimed 
crystals was acquired from (2), 
u = 94cm sec"!, 

being independent of their dimensions as shown in Fig. 1. One may see that 
the observed values of rimed crystals tend to increase with their dimensions. 
The disagreement between the values observed and calculated is perhaps due 
to the experimental formula (3’), in which the mass of rimed crystal is 
regarded as proportional to its area, being independent of its dimension.. 
More numerous observations for the rimed crystals are desirable. 


§3. Falling velocity of snow crystal of spatial type. 


In this case also, the difficulty in the estimation of # and a is not 
avoided, but the formula (2) may be approximately deformed as follows, 


_ / 2h(c—p) _ | 2ho (4) 
4d Gi nie Ce 


Assuming that the shape of spacial crystal is a circular 


cylinder as shown in Fig. 2, the volume is calculated, hh 
= git] 


volume = —- @h, 
4 Fig. 2. Side view of 


spatial snow crystal. 


where d is the diameter of the base plane, h the height. 
In Nakaya’s paper, the mass of the crystal of this type is described as fol- 


lowing formula, 
m = 0.0010 d2 


Thus the density is determined. 


c= 


TT 


36 C. Magono 


where v is the volume. Using the density, 
0.008 
u =,/ —-—— g= 49 cm sec}, 
Ri Cy, Tp g 

which indicates that the falling velocity of spatial type is also independent 
of its dimension. As seen in Fig. 1, the value calculated here is a little 
smaller than observed one. This disagreement may be due to the estimation 
C;: 1.2. Because the thickness of spatial crystal is comparable with its base 
dimension, Cy must be selected smaller than 1.2 probably. 


$4. Falling velocity of snow crystal of needle type. 


By means of the stroboscopic camera, it was observed that the snow 
crystal of needle type falls keeping its long axis horizontally. When the 
crystal has reached its terminal velocity u, its weight is balanced by the air 
drag which is perpendicular to the side plane of the needle, thus, 


ae Wao—p)g=Ccedl ein 2 , (6) 


where d is the diameter of needle, 7 the length, Cc the drag coefficient of 
the circular cylinder. 


— / rd(ao— : 
Thus, 4 = Yor oe g- (7) 
Calculating the masses for six samples of needle type in Nakaya’s book», it 
was found that their densities were about 0.15 g cm-3, being independent of 
their lengths. And plotting the diameters 
of the needles shown in the book to the d= 0.08 
length as seen in Fig. 3, following rela- 


tion was found, 3 
d=0.08#. S 
From the result of the experiment by the rf 
wind tunnel®, it is known that the drag = 
coefficient of circular cylinder set vertically oo 
to wind is about 0.7 in the range of Rey- 0 1 2 3 
nolds’ number 192~108, 5 =o, Therefore (i Kee NEEDLE IN 


(7) is rewritten as follows, 
u=14x102/7T cm sec-l. 


The falling velocity calculated thus, is plotted together with the observed 


values in Fig. 1. With the agreement of this order, we should be satisfied 
in this case. . 


* Dr. Houghton?) assumed d=0.1~0.075 if 
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§5. Falling velocity of graupel. 


By means of the stroboscopic camera, and from the Author’s observation 
of graupels fallen on a mire, it was ascertained that the graupel of cone type 
falls keeping its top upward, and as well known, its bottom is approximately 
hemispherical. So, it will be accepted that the drag coefficient of the sphere 
was adopted to calculate the falling velocity of the graupel. In the terminal 
State, the weight of the graupel is balanced by the air drag, 

TT TT 1 

MU pg t ees te ee eee 
where d is the diameter, og the density, and C,; the drag coefficient of the 
sphere. Thus we obtain 


“= Ae) g. @ 
The drag coefficient of sphere varies with Reynolds’ number, but the falling 
velocity of the graupel is yet unknown in this case. To simplify the calcul- 
ation, the value of drag coefficient: 0.5 was adopted according to Newton’s 
estimation. As for the value, it will 

be later described in detail. Thus the i‘ 
terminal velocity wis described by the | 
function of the diameter and density 

only. In Fig. 4, u is plotted by broken 
lines as the function of d with para- 
meter o. As shown by the dots in 
Fig. 4, the falling velocities observed 
by Nakaya, situate between the densi- 
ties 0.12 and 0.14 g cm-3. Considering 
that Nakaya gave 0.125 g cm-3 as the 
mean density of the graupels, it is 
accepted that the falling velocity of 
graupel is well described by (9), ex- 
cept the smaller range. In the range 

smaller than 2mm, the observed val- ei i 4 pb _ He = § 
ues are relatively small. This differ- Fig. 4 

ence may be due to the assumption 

that the drag coefficient is constantly 0.5 even in such a small Reynolds’ 
number 2x102. By the measurement of Dr. Gunn®, the drag coefficient of 
the water drop in these Reynolds’ number is about 0.7. It appeas that if 
the increase of drag coefficient with the decrease of Reynolds’ number is 
considered, the difference will decrease. On the other hand, the effect of 
scattering of the density should not be neglected. As shown in Fig. 5, it is 
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apparent that the larger the graupel, the larger the density. To: discuss the 
falling velocity in detail, it is necessary the simultaneous measurement of 
the density and the falling velocity of individual graupels. 


§6 Falling velocity of snow flake. 


In the previous report?’, the author described that the snow flake receives 
the resistance proportional not only to its horizontal cross-section, but to its 
volume. There, the falling velocity was derived from following equation, 


+ TR(ao—p)g=(Arr?t+ Brrt) 3 pe, (10) 


where + is the radius of snow flake of spherical type. The constants A, B 
adopted in the report, were inadequate. e 

It may be natural to consider that B is the function of the density of 
the snow flake, because the term B73 shows the resistance by a part of air 
which passes through the snow flake. 
In order to determine the function, 
the density of the graupel through 
which air can not pass, was con- 
sidered. The minimum density of 
graupel observed hitherto, -is about 
0.05 g cm-3 as shown in Fig. 5. There- 
fore, it is reasonable that B is propor- 
tional to (0.05—c) gcm~3, that is, B 
peak = (0.05—c)y. In the case of non- 
TAKAHASHI 
Rn permeable spherical body, the drag 
coefficient is estimated to be 0.5 as 
descrived above. (Reynolds’ number 
103~104). So A was considered to be 
DIAMETER OF GRAUPEL IN MM described by following formula, 


Fig. 5 A = 0.5—(0.05—o)x ’ 


DENSITY OF GRAUPEL IN G-CM™® 


where x and y are constants of no dimension. To determine those constants 
by the results of observations and formula (10), the simultaneous observations 
of the density, the radius, and the falling velocity are necessary. But no 
simultaneous observation has been made for those elements. During recent 
two seasons, the author measured the density of various snow flakes. Their 
volumes and diameters in falling state were determined by the “ photographic 
paper method”, and their masses were determined by the increase of the 
weight of the photographic paper by falling of the snow flakes. The results 
are plotted in Fig. 6. One may see there the tendency that the larger the 
diameter, the smaller the density. This tendency appears to originate from 
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the fault of the photographic paper method, 
in which the volumes of all snow flakes 
are calculated from the traces coloured 
on a photographic paper, neglecting that 
the large snow flakes do not always take 
a spherical form but often oblate ones. 
From Fig. 6, the mean density of rimed 
snow flakes for which many observations 


rimed 


non rimed had been made, was temporarily assumed 


DENSITY OF SNOW FLAKE IN G-cM™ 


to be 0.02 g cm-3. Besides the density, two 
mean falling velocities 195 cm . sec~! and 
208 cm sec~! of oblate and prolate type at 
the diameter 2cm and 3cm, were adopted 
to determine x and y. Introducing these values to (10), the constants x and 
y were calculated, 


0 0.5 1 15 2 25 
MEAN DIAMETER OF SNOW FLAKE IN CM 


Fig. 6 


core Vi 20. 


Thus the falling velocity of snow flakes was semi-theoretically acquired, 


(11) 


SSNPS ーー 7 
3 3 p  035—3a+ 1—20a)r * 


Assuming that tha densities of the non-rimed snow flake and wet snow 
flake were 0.011 and 0.025 g cm-3 respectively, the formula (10) describes their 
falling velocities as shown in Fig. 7. In the figure, one may see that these 
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calculated values of spherical type place half way of those prolate and oblate 
type in the wide range. Considering the density shown in Fig. 6, it may be 
accepted that the assumptions descrived above were not so unreasonable. 
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More simultaneous observations of the size, density, and falling velocity of 
snow flakes are desired to dicuss the falling velocity more exactly. 

In conclusion, the author wishes to express his best thanks to Prof. U. 
Nakaya whose observation occupies the essential part in this work, and to 
Mr. T. Maruo who gave the author useful advices in the aerodynamics. 
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Investigation of the Size Distribution of Precipitation 
Elements by the Photographic Paper Method 


By 


Choji MAGONO 


Abstract 


It is considered that the size distribution of solid precipitation is more 
important than that of rain drop to investigate the mechanism of precipi- 
tation. In this work, the relation between the diameter of various precipi- 
tation elements and the records coloured by them on the photographic paper, 
was mainly investigated. 


a) Rain drop. 
‘From the calibration of the water drops at terminal velocity, the follow- 
ing relation was obtained, 
d=10-2D?, 
where d= diameter of rain drop in mm, D = mean diameter of wet area by 
rain drop on the photographic paper in mm. 
b) Snow flake and snow crystal. 


The outline of the horizontal shape of snow flake or snow crystal in 
falling state is coloured and recorded on the photographic paper, so we obtain 
d=2-tyD, 
where d= mean diameter of snow flake or snow crystal assumed to be 
ellipsoidal in falling state, D = mean diameter of area coloured on the photo- 
graphic paper in mm. 


c) Very wet or partially melted snow flake. 


Very wet or partially melted snow flakes collapse and spread on the 
photographic paper when they fall on it. Their thicknesses are roughly 
uniform, so their diameters in falling state are calculated from their records 
on the paper and their uniform thickness 


d=(15)thsD?, 


where d= diameter, 7 = thickness of snow flake and D = mean diameter of 
coloured area in mm. 
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d) Graupel. 


The graupel melts to a hemispherical drop on the paper, so its diameter 
equivalent to that of a spherical water drop is calculated as follows, 


ds=2- ED 
where d, = equivalent diameter, D’ = diameter of dark brown area in mm. 
Some examples of the traces by various precipitations on the photo- 


graphic papers and their size distributions are shown in some photographs 
and figures. 


Introduction 


The size distribution of precipitation elements gives us an important 
method to search the mechanism of precipitation, and to use the result of a 
radar echo. While, except the rain drops no attempt has been made to 
observe the size distribution of other solid precipitation. It is not impossible 
to observe the mass distribution of solid precipitation by the “ filter paper 
method” like the rain drop, but in the case of solid precipitation, the simul- 
taneous observation of its sort and size is also desirable. Specially, in the 
sleety weather, the distinction between the rain drop and solid precipitation 
is necessary. 

Two years ago, the author noticed that snow flakes left their traces on 
a photographic paper by the daylight without any photographic operations. 
Since then he has made it a rule to utilize this phenomenon for the obser- 
vation of size distribution of precipitation elements. It is the purpose of this 
paper to show the relation between the size of various precipitation elements 
and their traces on the photographic paper. 


§1. Colouring of photographic paper by precipitation. 


It is well known that a photographic paper is coloured without develop- 
ment and fixing when it is exposed to the daylight. When precipitation 
elements fall on the paper and melt on it, the areas wet by them are coloured 
different from the other unwet areas as shown in photos. 3. By this 
phenomenon, the traces of rain drops or snow flakes. are recorded on the 
photographic paper. i 

Generally, the areas wet by precipitating elements are coloured to ini 
and the unwet areas to red brown. But for a time, the colours at both areas 
become violet brown. As the result, the distinction between them soon 
vanishes. Fortunately, on a sort of photographic papsr* in our country, 
unwet areas are coloured blue, and wet areas red brown, when exposed to 


* “Yae FS”, “Somei FS”, “ Yoshino No. 4”, “ Hinode No. 2”, made in Konishiroku’ 
Photo. Ind., Ltd, 
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the direct daylight or a light of relatively short wave length, as shown in PI. J, 
Photos. 1~3. In addition to the clear difference between the two colours, these 
colours are stable if the paper is protected from the direct rays of the sun. 

In this work, photographic papers “Hinode No. 2” were mainly used. 
The direct daylight, brilliant electric lamp of magic lantern, and mercury 
lamp etc. are suitable for the light source of this purpose. The exposing 
time ranges from several seconds to ten minutes, according to the intensity 
of the light used. After the paper was coloured properly, it must be natural- 
ly dried and preserved at a place protected from the moisture. As this 
method needs no photographic operations, it is very suitable for the field 
work. Even an aged or somewhat fogged photographic paper can be used 
for this purpose. 


§2. Observation of rain drop. 


It is necessary to calibrate previously the relation between the size of 
the trace coloured by a rain drop and its diameter. 


a) Calibration. 


Because the size of the trace by the water drop varies with its falling 
velocity, the calibration must be made with a water drop of terminal falling 
velocity. According to Laws report,” the height of falling requisite to attain 
to terminal velocity was selected as in the following table. 


Table 
Diameter mm Height m 
0.3~0.6 3 
0.8~2 6 
3 ~6 8.6 


The drops’of diameter 3~6mm were produced by dropping water from 
the end of the glass tube as shown in Fig. 1-A. Their diameters were 
determined by the increase of weight of a photographic paper due to the 
fall of the drops. Such a large drop leaves a brown trace and a dark brown 
trace in the brown one as shown in Photo. 1. It is supposed that the drop 
spreads in the brown area when it falls on the photographic paper, then 
shrinks to the dark brown area. In this work, the mean diameter of the 
brown dark area was adopted as the size of the trace. It is difficult to 
receive the drop of diameter smaller than 3mm with the paper by the 
“dropping method,” so the spraying method was used as shown in Fig. 1-B. 
If the water pressure in the nozzle is stable, the distribution of the drop size 
is invariant. The diameter of those small drops was calculated from the 
_mean weight of many drops of almost equal size which had fallen on the 
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paper coated 
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Fig. 1 


photographic paper. The trace diameter was measured from the mean of 
traces on the photographic paper with a magnifying glass. These traces are 
shown in PI. IJ, Photos. 4,5. The drop size smaller than 0.6 mm was determined 
by the “ microscopic method” in which the side section of the drop which had 
fallen on a paper was measured with a microscope. For this pupose, a paper 
coated by paraffine was more suitable than the photographic paper, because 
the side section of a drop on the photographic paper is too much complicated 

to calculate the volume from its side section as shown in Fig. 1-C. The size 
of their traces was observed by the other photographic paper on which other 
drops are of equal size to that which had fallen on the coated paper. See 
Photo. 7. 


b) Relation. between the size of trace and rain drop. 


Using the apparatuses descrived above, the relation between the trace 
‘diameter D (mm) and the water drop diameter d (mm) was acquired, and 
the result is presented by a full line in Fig. 2. The broken line in the figure 
is the calibration curve for the “ filter paper method” given by Takahashi, 
The curve coincides approximately with that for the -photographic paper 
method, except the range of diameter smaller than 3 mm. This fact may be 
explained as follows. The thickness of the water drop shown by the brown 
‘area on a photographic paper coincides accidentally with that of filter paper 
saturated with water, while in small drop range the extension-of water by 
‘falling on the photographic paper, does not reach to that in the filter paper, 
‘because the momentum of falling is small.- The size of the brown trace-over 
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which a drop has spread by falling is not considered to be related with the 
sort of the photopaper except in the case of very small drops. The trace 
sizes on the photographic papers of different sorts were. practically equal 
each other if only smooth photographic papers were used. Considering the 
stability, the clearness of colour, and from the economical point, photographic 
papers “ Hinode No. 2” (Gaslight) of thin type were suitable for this pupose. 

Changing the ordinate of Fig. 2 to logarithmic scale, the result of the 
calibration is shown by a straight line which passes through the origin as 
seen in Fig. 3, and is formulated as follows, 


log (10 D)=1.5 log (10 d) 
d=10-4DI=046 Dt, (1) 


where d is the diameter of a rain drop, D mean diameter of a trace by a 
water drop. This relation corresponds to the relation between the area and 
the volume. In other words, it is shown that a water drop spreads in a 
moment at a uniform thickness on a photographic paper. The uniform 
thickness #-(mm) is calculated from the following equation, 

“pPh= volume = 3, - , eos 


oo i nh - gat) 43 
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From (1) and (2), h is given by 0.067 mm. The error of the diameter by 
this method is supposed to be +4%, considering that the difference between 
the total mass calculated by this method and the one measured by weighing 
is +10%. 


c) Observation of rain drop distribution. 


Photos. 6, 7 in Plates II, III, 
show the traces of rain drops 
on photographic papers obtain- 
ed at Shiozawa, Niigata Pre- 
fecture, in Jan. 1953. Figs. 4 
and 5 are each size distribu- 
tions from the photographic 
papers. Photo. 7 in Pl. HI and 
Fig. 4 show the narrow size 
distribution of small rain drops 
at the end of rain fall process. 
In Fig. 5, one sees two maxi- 
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medium drop size, as Taka- 
hashi® pointed out. The capital letter N in these figures shows the number 
of drops observed. 
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Fig. 5 
§4. Observation of solid precipitation. 
a) Snow flake or snow crystal. 


It may be considered that a trace by a snow flake corresponds to the 
horizontal cross-section of the snow flake in falling state, because the ordinary 
snow flake falls on the photographic paper without spreading, even if it 
collapses somewhat longitudinarily. Assuming: that the snow flake takes a 
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rotating ellipsoidal form in falling state, the volume of the snow flake is 
measured as follows, 


volume = =s abv ab. (3) 


where a is the major axis of the trace, b the minor axis, 
and ab is regarded as the diameter of vertical axis of 


snow flake in falling state. See Fig. 6. Yab is also con- 


sidered to be the mean diameter of the snow flake, . 
Fig. 6. Trace of 


a snow flake 


d=(abVab)*=V ab . (4) 


By measuring the increase of the photographic paper by the fall of snow 
flakes, the mean density of the snow flakes which has fallen on the photo- 
graphic paper is calculated. The snow crystal was also treated as the small 
snow flake. The weight of the photographic papers used, was about four 
grams, and the increase of the weight was of the order of a few hundred 
milligrams. 

In PI. I, Photo. 3, an example of the traces of the snow flakes obtained at 
Shiozawa in Jan. 1952 is shown, and the size distribution is shown in Fig. 7. 
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It is supposed that the larger maximum in the volume distribution shows 
the distribution of the snow flakes, and the smaller one shows the distribu- 
Hon of the single crystals and the fragments. The followings are the charac- 
ters of snow flake distribution, that the volume distribution is entirely 
different from the number distribution comparing the case of rain drop, and 
that in the range of large diameter, the distribution is very irregular, 
although the observed snow flakes are numerous enough. For example, the 
volume distribution 22.4% in the range of diameter 17~18mm was occupied 
by only one snow flake. The total volume .and the total mass of the snow 
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flakes’ -on -the- photographic were 18cm’ and 280 mg -respectively; therefore 
the mean density was 0.016 g cm-3. 

In the case of those snow flakes which are scattered when they fall on 
the photographic paper, it is difficult to determine their volumes. But the 
mass distribution is obtained by weighing the photographic paper, for the 
fragments produced from one snow flake are assembled usually to some 
extent. : 


b) Very wet or partially melted snow flake. 


A very wet snow flake spreads horizontally to an uniform thickness 
when it falls on a photographic paper. But it does not shrink by surface 


tension like a rain drop, because it contains many ice fragments unmelted. - 


The thickness of the wet snow flake on the photographic paper varies with 
the grade of its water content in the range from 0.067 mm to about 1mm. 
Fortunately the thickness is roughly common to all snow flakes which have 
fallen at the same time, because their wetness is similar to one another. 
Measuring their uniform thickness, their volumes and mean diameters are 
determined as follows, 


volume == eB: Dh= ae a3 


— Vi5h D'=115h*D*. (5) 


When the snow flakes are almost melted, therefore, their densities are 
considered to be about lg cm-3, their common thickness # is calculated by 
weighing the total mass on the paper, 


h=M/2 D, (6) 


where M is the total mass of the snow flakes. 

Fig. 8 is the result of wet snow flakes obtained at Sapporo in April 1951, 
when the thickness was uniformly about 1mm. By comparing the total 
volume thus calculated with the total mass by weighing, their mean density 
0.7 g cm-~3 was obtained. 


©) Graupel. 


' Some time after when a graupel falls on a photographic paper, it is 


melted into a water drop of hemispherical form, as shown in Fig. 9, In this 
case, the dark brown area corresponds to the base of the hemisphere. Thus 


the equivalent diameter of the graupel to the orate of spherical water 


drop is obtained as follows, 
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where d,, is the equivalent diameter of the graupel, D’ the mean diameter 
of the area of dark brown colour. The volume of a graupel can not be 
determined by this method. The error in the equivalent diameter of graupel 
is supposed to be +3%, comparing with the absolute weighing. 

An example of the graupel which fell at Yokohama in Jan. 1953 is shown 
in Photo. 9, and the distribution is presented in Fig. 10. In the figure, one 
may see that the size distribution of the graupel is similar to that of the 


rain drop. This suggests that the graupel directly becomes the rain drop 


when it melts. 
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§4, Rain with solid precipitation. 

a) Rain drop and very wet snow flake. 

When rain drops fall with very wet snow flakes, this photographic paper 
method is very convenient to distinguish them. As seen in Pl. III, Photo. 8 
which was obtained at Shiozawa in Jan. 1953, almost melted snow flakes are 
easily distinguished from rain drops, for the forms of the traces by the wet 
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snow flakes are complicated, and some splashes by unmelted ice pieces are seen 
in the traces. In such a case, the equivalent diameter of almost melted snow 
flakes is calculated by following considerations. By reducing the sum of the 
volumes (masses) of rain drops from the total mass of the rain drops and 
the snow flakes, the sum of the masses of the snowflakes is obtained. Using 
(5) and (6), the thickness and the equivalent diameter are calculated. 

The size distribution shown in Fig. 11 was thus obtained. In the figure, 
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one sees that the size distribution of very wet snow flakes deviates to the 
larger range, in other words, small snow particles, such as single snow 
crystals or fragments had melted earlier than the snow flakes. This is a 
type of the sleety weather. 


b) Rain and wet graupel. 


In Dec. 10, 1952, many ice particles were observed on a photographic 
paper at Yokohama, being mixed with rain drops, though these ice particles 
were not noticed until they had fallen on the photographic paper. The 
particles remained for some time unmelted, and the traces were ‘coloured 
more powerfully than that of rain drops. In order to make sure, the arrows 
were marked to their traces as shown in Photo. 10. The author thinks that 
the particles were wet graupels, for they were not transparent. And if the 
particles originate from refrozen rain drops, that is, sleets they must be 
smaller than the rain drops. For the freeze must begin from the smaller 
‘rain drops. But their volumes (calculated from semispherical form) were 
‘rather larger than that of rain drops. 

In Fig. 12, the size distribution of the rain drops and graupels which were 
obtained from Pl. IV, Photo.10, are shown. Each diameter was-calculated by 
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Photo. 7 (Fig. 4): x1 
Rain drops. 
A part of Cabinet size, expos. time: 2 min. 
13h20m, Jan. 25, 1953, at Shiozawa. 


Photo. 8 (fig. 11): x1 
Rain drops and wet snow flakes. 
A part of Cabinet size, expos. time: 10 sec., 
13h10m Jan. 25, 1953, at Shiozawa. 


Photo. 9 (Fig. 10): 


Graupels. 
A part of Cabinet size, expos. time: “30 sec. 
23h15m Jan. 7, 1953, at Yokohama. 


Photo. 10. (Fig. 12): x1 ; 
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marked by arrows. , 

NE of Chl Snes 30 sec. 
OO De EE ee 


Investigation of the Size Distribution of Precipitation Elements 51 


each calibration method. In this case also, the graupels distribute in the 


range of diameters larger than those of rain drops. 
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Conclusion 


A simple method by the photographic paper which is applicable to all 
precipitation elements, was proposed. But the method has yet following 
faults. In a heavy precipitation, the traces of large rain drops mix mutually. 
The volumes of those snow flakes which scatter on the photographic paper 
can not be measured. These faults are common to all other methods in 
which the precipitation elements are received by a paper, but they shall be 
removed in the near future. In a rather temperate climate, the author 
thinks, this photographic method is most convenient. 

Recently the author has known that Dr. Smith had used already a 
sensitive paper for the approximate measurement of the rain drop size on 
an aircraft, but his detailed procedures have not been reported yet. 

In conclusion, the author wishes to express his best thank to Misters. 
B. Arai, T. Watanabe and M. Matsumoto for their collaborations of this 
work. The author is also indebted to Mr. M. Shida who gave him the 
facilities for the observation in this work. 
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Abstract 


The well known divergences due to the energy conservations in the 
intermediate states are considered, from the point of view of the graphical 
structure and of the field reactions. 


§1. Introduction 


In contrast to the vitraviolet divergence, the low frequency divergences 
due to the coincidence of poles in transition matrices or in transition probali- 
ties in the quantum field theory are generally regarded as a defect of the 
perturbation theory. These difficulties should disappear in principle. How- 
ever, they stand sometimes in the way of computational works with the 
Feynman-Dyson’s covariant S-matrix formalism, especially in the proc2ss 
involving mesons. Therefore, the special discussion of these divergences is 
also a matter of practical importance. 

In quantum electrodynamics, the low frequence divergences are conveni- 
ently classified into three cathegories; the first of them appears from the 
single pole of matrix elements corresponding to the energy conservation in 
the intermediate states. As is well known, the divergences of this type are 
of quite formal nature and they raise no difficulties, as they are removed by 
well known procedure of the amplitude-renormalization of the wave functions, 
or by other suitable limiting processes, such as the principal value evalua- 
tion, etc. Therefore, they are out of our consideration. The second type 
divergence originates from the matrix element which is finite certainly but 
not uniformly. A typical example of this divergence is so called “ infra-red 
catastroph”, that is, the matrix elements of bremsstrahlung are certaintly 
finite, but their total transition probability diverges at the low frequency 
limit. The divergence of the third kind comes from the coincidence of poles 
in a transition matrix. As is well known, the coincidence of the displaced 


poles corresponds to the excitation of the real process or to other competing , 


processes. Therefore, the divergence of this type is regarded as a “ resonance 


*) The content of the present work was published in Japanese more one or two 
years ago, separated to several papers. So, we summarize here our. main considerations. 
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catastrophe” and they will disappear whenever the effects of field reactions, 
i.g., the effect of radiation damping, etc., are properly taken into account. 


In the meson theory, it is evi- 
dent that no infra-red catastrophe 
can appear on account of the non- 
vanishing meson mass. The reason 
why this is so is readily seen if we 
note the fact that the infra-red 
catastrophe originates from the 
vanishing energy denominators in 
the perturbation theory. As a simp- 
lest example we consider the process 
of bremsstrahlung in Fig.1. As the 
real emission of a quantum & on 
the transition from free particle 
state p into a free state p—k is 
prohibited by the conservation law, 
the matrix element of this process 
is certainly finite, except only for 
one case k=0, where the energy 
denominator of the intermediate 
propagator vanishes. The element 
is finite in almost everywhere, 
however, it is not uniform. But in 
the corresponding meson process, 
k=0 is excluded by non-vanishing 
meson mass. This is the reason 
why no infra-red catastrophe ap- 
pears in the emission process of 
meson. 

But in the case of the corres- 
ponding process of meson-decay 
where the meson pf in intermediate 
state decays into two lighter parti- 
cles, a free particle state is allowed 
as the intermediate state. There- 
fore, new divergence corresponding 
to “infra-red catastrophe” re- ap- 
‘pears in this decay-process. Let us 
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Fig. 1. V: external field of force 
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Fig. 4 


now consider this “infra-red-like catastrophe” divergence. For this purpose, 
it is convenient to introduce a Feynman-diagram of transition probability. 
This diagram is constructed from two ordinary Feynman-diagrams G, and 
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G», as follows. We define here the conjugate diagram G* of G, which is 
obtained by reversals of all direction of the propagating particles in G. Put 
G* at the right edge of G; and join each lines smoothly as illustrated in 
Fig. 2, then we obtain the Feynman-diagram of the transition probabilities 
constructed from G; and G2. In the same way, we have Fig. 3 as a Feynman- 
diagram of, a radiative correction term for the transition probability of the 
elastic scattering by an external field of force. Here it is notable that Fig. 
2 has the same structure as that of Fig. 3 with different position of the 
“line of observation” (+00). A part of the transition probability, which 
results from a cross term of the matrix elements H, and Hs corresponding 
to the graphs G; and Gy» respectively, 


W= “7 <H, H;> dd) 


is readily obtained from these diagrams according to the usual rules of 
Feynman, provided that the “ internal” lines which cross the line (+00) in 
Fig. 2 must be regarded as the out-going free waves in the final state, that 
is to say, the observation- (or projection-) operators to the final states are 
hidden on this line, 

The problem of the infra-red catastrophe has been discussed by many 
authors,” and clarified that the infra-red divergence appeared in the brems- 
strahlung is cancelled out by that of the radiative correction of the elastic 
scattering process. Recently, T. Kinoshita?) analyzed this problem in more 
detail and demonstrated that cancellations are accomplished term by term 
between the parts of transition probabilities that have the Feynman-diagram 
of same structure. Extending his result, it is shown that the same procedure 
can be successfully applied to the case of the infra-red-like catastrophe diver- 
gence in the meson-decay process and finite rate of decrease of initial state 
probality is obtained in each order approximation of perturbation. This 
method gives, however, no finite results for each one of the component 
process. This situation often limits the practical applicability of this method, 
inspite of its theoretical interest. In order to obtain the finite result for 
each component process, we have to take into account the effects of the field 
reactions, such as the effects of radiation damping. The effect of radiation 
damping is introduced here entirely within the frame-work of Feynman- 
Dyson’s covariant theory of S-matrix. A covariant generalizatien of Heitler’s 


theory of radiation damping has been proposed by T. Miyajima and N. 


Fukuda? So far as the effect of radiation damping is concerned,. our 


1) F. Bloch and A. Nordsieck, Phys. Rev. 52 (1937), 1045 (L) 

W. Pauli and M. Fierz, Nuovo Cimento 15 (1938), 167. 

W. Braunbek and E. Weinmann, Z.S. f. phys. 110 (1938), 360. 
2) T. Kinoshita, Prog. Theor. Phys. 5 (1950), 1045 (L). 
3) N. Fukuda and T. Miyajima, Prog. Theor. Phys. 5 (1950), 849. 
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method gives, of course, equivalent results with that of Miyajima and Fukuda. 
The difference consists in the way of taking into account the virtual processes. 
As our method is developed completely in conformity with the Dyson’s S- 
matrix theory”, it is convenient for the renormalization, for S-matrix theore- 
tical characterization and for application to the practical problems. 


§2. Infra-red-like catastrophe aspects.>.© 


In this section, let us show that the infra-red-like catastrophe divergence 
which appeared in the meson-decay process is removed by the same procedure 
as used in the theory of infra-red catastrophe. In order to clarify the cor- 
respondence to the ordinary infra-red catastrophe, we consider the meson- 
analogue of bremsstrahlung, i.e. the process in which a heavy meson (mass 
M) decays into two lighter mesons (mass m and nm) in the course of deflec- 
tion by an external field of force V. (See Fig. 1) In contrast to the ultra- 
violet divergence, the type of the relevant particles is of minor importance. 
Therefore, we may choose the simplest model for each particle without any 
loss of generality. As an illustrating example, let us consider a part of the 
transition probability, the Feynman-diagram which is illustrated in Fig. 4. 
Using following the simplest model ; 


Ux) : charged heavy scalar field of mass M 

Cx) : charged light scalar field of mass m (2.1) 
P(x) : neutral scalar field of mass m/’ 

V(x) : external field 


and their mutual interactions ; 


H(x) =g(Up+ yl) (2.2) 
K(x) =U UV (2.3) 
KD=¥ YYV (2.4) 
the transition probability Wa of the decay process “A”™ is calculated as 
goa ale (ae EO FE P 
4(27)8 Es sg Cp Ebr-k CEs—E $7] Ces-n—ep-z] 


x 6(Ep—ep/-k—Ex) (25) 


where (£4, €p, €p) =/ P2+ (M2, m2, 2) respectively. The integand has two 
poles. As is easily seen, one of them corresponds to the energy conservation 


4) F.J. Dyson, Phys. Rey. 75 (1949), 1736. 

5) D. Ito, Prog. Theor. Phys. 6 (1951), 1020 (L) and 1022 (L). 
6) J. Edden, Proc. Roy. Soc. 210 (1952), 388. 

*) The process “A” in Fig. 4 is the same with Fig. 2. 
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in the intermediate state B, in Fig. 4, and the other corresponds to the energy 
conservation in the intermediate state C. In the vicinity of these poles, Wa 
diverges logarithmically and takes the following asymptotic forms. 


1 dK dP’ T A sem 2 8(Ep—ep-kr—€Ex) 
W4(B)= (2.6) 
Ae ‘cb > Es 2Ep, ep-kEp Cb-p— Ob-k Cp-k+ Er— Ep 
in this case (€k~E yp —ep-pz), 
_ gt 1 (_ dK aP | VCP=P)P (Ep SEED 
as 4(27r)8 TA serpin es E}—E, € p/-k—Cp-k Son, 


in this case (€p/-k=€p-p) 

In analogy with the theory of ordinary infra-red catastrophe, let us 
examine the transition probability Wg of the elastic scattering process “B” 
in Fig. 4*™ In this case, we have to assume the line B as the line of the 
observation. Ws is calculated as follows; 


Ws= 


£ 7. 上 Sil dP’ | VC(P—P') #2 
2m EpJe,-u ek J Ep Ces -4—(Ep—e, Wey CE —Ek] 
x 6(Ep—Ep) (2.8) 
Two poles of the integrand correspond to the conservation of energy in the 


intermediate states A and B respectively, and the asymptotic forms at these 
points are 


ET DE | dKdp “FWP-P)E 8(E,—Ep’) 
Wr(A Da | Dp D i 
eS 4(2ar)8 Ep) 2Ep ep re, Cp-k—Cp-k Ep-k+ Er— Eg 22) 
in this case (€g*-Ep—ep-~) 
eee dKdP |VP—P) 2 8(Ey—Ep ) 
Wa Lda Es 
BC) es, Ey) 2Ep ep-k Ek br Cp-pt+er—Ep (2.10) 


in this case (ex +Ep—ep-, ) 


Further, we have to take into acount another process whose Feynman-diagram 
is obtained from Fig. 4, assuming the line C as the observation line. This 
is a process of meson-decay accompanied by the double scattering by the 
external ‘field, which has no analogue in the corresponding photon process. 
The transition probability We of this process ““C” turns out to be 


Wer ("eK (ae | VP—P') 2 
(Sts 4(2 8 E 2 —E?2 2 2 
(27)8 Ep A: Ee J ep-k LEp—Ep]lep-.—e-x] 
ko 
x 6(Ep—ep-7—E ) (2.11) 
and the asymptotic forms at A and B are 


**) The process “B” in Fig. 4 is the same with Fig. 3.. 


i = es re 


ee eee 
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Wc(A)= £ 8 : I amir SAE eae 2 S(Ep—ep-4 eg) 
ZT Ey) 2es-pep—p Ey Ey—E} Cp! -p—Cb-p 
in this case (@p/-p2=€p-, ) Cay 


WAB) = —& a GK GPT, CP nh) BSCEsereps a7 66) 09,43 


4(27)8 Ep) CIE EL EC EC Ep/—p-p—€p 
in this case (é€z=Ep, —e€p-g) 
Comparing these asymptotic forms of divergent integrals, we can readily 
find out a kind of “reciprocity relation ” holding between them, 
Wa(B)=—W2(A) 
Wa(C) =—Wc(B) 
Wc (A)=—WaA(C) (2.14) 


From the view-point of the asymptotic forms of the divergences, we may 
summarize the above results as follows 


Wa~ WaCB)+ WaCC) 
Wa~ W2(A) +WRC) (2.15) 
We~ Wc(A)+ We (B) 


On account of these asymptotic forms and the above reciprocity relations, 
the infra-red-like catastrophe divergences are cancelled out from the total 
rate of decrease of the initial state-probability. 

In this way, the divergence appeared in the process of our meson-decay 
is cancelled out by simultaneous consideration of the processes of the elastic 
scattering and of the decay with double scattering. This is nothing but the 
way, by which the ordinary infra-red catastrophe has been removed. There- 
fore, the divergences due to the energy conservation in intermediate states 
of the meson process are of the same nature as that of the “infra-red 
catastrophe” and they are removed by natural generalization of the method 
used in the case of quantum electrodynamics. 

Moreover, the above consideration reveals us that the reason why the 
cancellations of the infra-red-like divergences are possible, consists in the 
validity of the reciprocity relations. Of course, these relations also have 
played essential roles in the case of the ordinary infra-red catastrophe, but 
we could not notice them on account of the over-simplifications of the mecha- 
nism due to the simpler nature of photons. Accordingly, our method will be 
extended to more general cases, if we are able to give a general proof of 
the validity of the reciprocity relations. In order to give it, let us consider 
a Feynman-diagram of the transition probability of more general form, ‘as 
illustrated in. Fig. 5, where each one of G;; Gz and G represents arbitrary 
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graph. Referring to this graph, the reciprocity relation to be proved is 
formulated as follows; The asymptotic form of divergence W4(B) due to 
the energy conservation in intermediate state B of the process “A” has just 
the same form and opposite sign as the asymptotic form Wg(A) of similar 
divergence due to the conservation in A of the process “B”. Assuming 
again all relevant particles are scalar, we obtain the following asymptotic 
forms after some elementary calculations 


a GP; 2 ( dq; fi, ECP) 3 dC CT 3) 
py wt aa EP) 5d 
Wa®) fi) E;CPi) 35-1) €e;(qgj) > eq)—D ExT) " 
x 8(Einitial — = e;(qj)) (2.16) 
fz 
in this case (> €;+=>) Ei) ご 
& a aq; FP; Ei gd eh) 3) 
AN lee _ ( x 
Be a Ei(P;) i=1) €;(qj) > Fi) ES ICT) 
x 6 Ejnitial — DD (2.17) 


in this case .(3} => e) 
These results show that the reciprocity relation . 

Wa(B)=— Wa (A) (2.18) 
holds. for the most general cases. The possibility of the cancellation of the 
divergences due to the energy conservations in intermediate states by means 
of simultaneous consideration of related processes is an immediate conse- 
quence of this reciprocity relation. Accordingly, the method once applied to 
the cancellation of the infra-red catastrophe is now able to be extended to 
more general cases of such low frequency divergences as arising in the 
meson theory. 


(-«) A B (-%) 
Fig. 5. ECPD=V PHM? s(a,=" + 
§3. Consideration of field reaction 
In the previous section, it is clarified that the method used in the 


theory of infra-red catastrophe can also be applied to cancellations of the 
low frequency divergences appearing in the meson-decay processes. However, 
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we could have only a finite rate of decrease of the initial state-amplitude. 
The transition probability of each component process which transfers the 
system in its initial state into each final state is still divergent, therefore, 
we could have no information about these component processes. This. situa- 
tion limits considerably the applicability of this method, inspite of their 
theoretical interest. In order to obtain a finite result for each component 
process, the effects of field reactions must be taken into account. 
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(1) general survey 
Let us begin with a simple example; the double decay process of mesons, 
in which a heavy meson (mass M)-decays into two lighter mesons (masses 
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m and m’) and then one of them (say, mass m,) decays into two lighter 
mesons (mass » and pw’). The Feynman graph of this process is illustrated 
in Fig. 6. According to the perturbation theory, the lowest order approxi- 
mation of the total transition probability diverges at the point which cor- 
responds to the free particle state of the intermediate particle m. As men- 
sioned above, the divergence of this type is of the same nature as that of 
the infra-red catastrophe in quantum electrodynamics and will be removed 
by the same method. 

Apart from the divergence-problem, however, the description by the 
lowest order approximation of perturbation is expected to be qualitatively 
very incomplete. Because, unstable mesons are fated to decay sooner or 
later, and they will certainly decay into lighter mesons in a very long lapse 
of time. And whenever the rate of decay of the intermediate meson m 
is greater than that of the creation, the rate of the total process will be 
determined by the decay probability of the mother meson M alone, that is, 
the S-matrix element of the total process M>m+m’'—>m'+p+p’' will be in 
the limiting case independent of the detailed forms of the interaction operator 
describing the latter process m—>yu+ yp’. It is probably impossible to describe 
such situations by the lowest order approximation of perturbation alone. 
The above situation is clearly a result of very long duration of time after 
the creation of meson m in comparison with their life. The transition proba- 
bility W~|<pp! | H,|m><mm'|H,|M> |t calculated to the lowest order 
approximation has its meaning only if the duration of time is very small in 
comparison with the life time. Nevertheless, if we tend ¢ to infinity in order 
to know about the situation in infinite future, W increases beyond all limits.” 
This is nothing but the space-time interpretation of the infra-red-like cata- 
strophe divergence, as the linear time dependence is a result of the energy 
conservation in the intermediate states. Therefore, the divergence will dis- 
appear in such theory as is able to describe correctly the behavior of transi- 
tion probability of the total process. 

From above considerations, it seems natural to remove the divergence 
by considering the effect of damping in the intermediate states. In view of 
the characteristic features of damping, it is expected that the situation will 

‘probably be described by the following natural interpolation, at least as an 
approximation. 
W~|<mm' | H,| M> |? C(t» life) 
>Wel<M Imm > Px 
W~|<pp’|H2|\m><mm'|H,|M>\t (é«life) 
xX (1—e-|< m| Hal mm’ > Pt (3.1) 


This is not an unfounded expectation. When the proper frequency of the 


7) B.A. Lippmann and J. Schwinger, Phys. Rev. 79 (1950), 469. 
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meson is larger than the uncertainty of energy due to the finite life time 
2 : , ae 
rant that is, when the particles in intermediate state can be regarded as if 


it were in its free particle- state, the whole process may be a series of the 
natural radioactivity. In the series A, B, C, of natural radioactivity, the 
probabilities of finding the atomic nuclei A, B, C, respectively, 74, mg and 
mc are solutions of the following equation 


Na=—I'na, he=—ynptl' na, c=+y Rp 
with initial condition 


na(0)=1, xp(0)=nc(0)=0 


where /’, and y represent the transition probabilities of the process A っ 5B, 
and B-C, respectively. The rate of the total process W = 7c is calculated as 
W=I'(—e-") (3.2) 
This is just the same form as (3.1). Therefore, the correct theory has to 
have a solution which reduces to (3.1) when the intermediate state could be 
regarded as a state of free particle. Such solution is really obtained if we 
take into account the radiation damping in intermediate state in accordance 
with the theory of Weisskopf and Wigner. 


(2) damping effect in intermediate state 

Suggested by the above rough consideration, let us consider the correc- 
tions of the internal line by succession of the self-energy loops consisting of 
the pair of lighter mesons mw and w’, as illustrated in Fig. 7. Assuming all 
the particles are scalar and, their interactions are 


Hy=2.*6¢'+conj. H,=9,$* aa'+conj. (3.3) 


where W, ¢, ¢’, a, and a’ are fields of meson M,m,m’',u and yp’ respectively. 
We obtain the following expression as total contribution to S-matrix from all 
graphs having the forms given in Fig. 7. 


= 181& 1 : (3.4) 
(2m )3 CLA EC ob tm —ilet+ FH) 


(dk"’) (3.5) 


= & 
where P(p)= Cm), [R24 we—te L(p—k" 2+ p27] 


represents the self energy operator of the internal line with propagation 
vector p. The transition probability W is calculated from above S-matrix as 


(21 dP dK 1 = 
eter | OC rr Ca Te a 


8) V. Weisskopf and E. Wigner, Z.S. £. Phys. €3 (1930), 54 and 65 (1930), 18. 
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In the decay process, J’(p) is, in general, complex and given by 


aT (p) = 3m) + ify YP 8.2) 
2 OC(h—B!)2+ pL) . OCR? + w?) ) 
aaa 2 i 
= ti |e) (ort pepe 
+ im? 8(( p—BP-+ wl) BCRI2+ | (3.8) 


where the real part 62(p) represents the self-energy of the intermediate 
particle containing a correction of mass, an effect of level-shift etc., and was 
a principal matter of concern in the theory of the renormalization. The 
imaginary part ho ais contains in this case a non-vanishing part 
adK/(27)3 
8 fp Ex" Po— Ex") 
As is readily seen this is nothing but the total transition probability of the 
second decay process m—>y+p’. Putting €p=V P+ (m?+8m2(p)) we have 


= (£1 82)? | aP dps dK Ce + €p-4—Do) 
(27)5 J 16E Ep—pyoerep-2 Cp— 25)2+ (Doy(b)/2Y 


As y(p) does not vanish, we obtain always convergent W. The functional 
forms of 6722(p) and y(p) are usually very complicated and it is not a simple 
matter to evaluate the integral rigoroutly. It is not our present purpose to 
calculate the transition probability exactly, but to verify the fact that the 
result is divergence-free and is insensitive to the detailed nature of the 
interaction describing the second decay process, whenever the free particle 
state is dominant in the intermediate state and the latter decay is much 
faster than the first. Therefore, we calculate it here approximately in the 
case of “sharp resonance” where y(p) is small enough to replace the argu- 
ment p of dm(p) and y(p) by its resonance value. Then, the evaluation of 
the integral is elementary and reduces to 


W=2 i| CaP)/(2T3 | dK/(27)3 
= BE, Ep—ep)ep an Sep, (ep—ep) 


x 5(ep—En—Ep-z) WP) (3.11) 


yp) = 2m gi | Ce (3.9) 


W (3.10) 


On account of the definition of y(p), this reduces further to 


CaP)/(2T)3 


BE CFD) 2 


as 2rgt | 


This is nothing but the transition probability of the decay process of the 
first step M>m+m’. . Thus, under the condition described above, the transi- 
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tion ‘probability of the total process is determined by the first decay process 
alone, and does not depend upon the detailed natures of the second decay as 
expected before. 


(3) time- dependence of transition probability 


The above results are based upon the S-matrix theory and are able to 
predict only the situation in infinite future. The next task is to examine 
the time dependence of the transition probability. 

The transformation function Ufo, oo] of this process is written as” 


UEo.o5 |= 3 a dx! dx, dx! dn dxn! A(x) A(x—x!) L(x! —x) x 
n=lJ eo 


XxX Ai—L!) LC — 2): X LC 1— 4 ) ACtn—%n!') Br!) (3.13) 
where 


A(x) = 1 C4) Wx) LC%) =—(go/2)? Dp(x) Dix) 
BCX) = g2 a(x) a(x) CE =< Pa), dX) 5. Gil) 
A(x) = 4 Ap(x) 4 Dr(*—#/) = < P(a'(®), a'(x')) >o 


or in momentum space, it turns out to be 


Ula, 013 | dp dps dpa dpm ASD) ACP) LG,Cb, 1) AC.) x 


x tLe. Pn=19 Pad A( pn) Be,(bn) (3.15) 
Cs 1 cts idx sf Pp OE 4 — ~ipxtri 
A any | A@ ede, LD = gt | LD eitrtirdr dy 3.16) 
Be | ee e-ibtdx Aes 
70 2? Joy ip? + m>—ie) 


As the contribution from each one of the loops L7,(f,q) can be written as 


12,9) = eit-orasl 1G) e-i9¥ dy = 82,(p—q) Le--O(q) (3.17) 
0 


To 


(274 


we have 
Ulo, a0] = 3a dp, dpn SG, (bp—b) “826 Pn-1—bn ) Az CD) 4(b) x 


x LG-0(p;) Ap) LE (pn) A( pn) BE Chm) (3.18) 
where 


NO ' WIS Neal rere gear 
dx ~Ub-D« = §(p—q) a ebo- Wt’ dt (3.19) 
(27r }4 | 2 21 to 
reduces to an ordinary 8-functions when the time distance t—f)= T is large 
enough. As we are interested in the effect of the radiation damping, it is 


(hp—g) = 


9) D. Ito, H. Tanaka and M, Yamazaki, Prog. Theor. Phys. 7-(1952), 128 (L). 
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of no use to consider the behavior of the transformation function U[e, oo] 
in such a region where the temporal distance T is comparable with the life 
time Consequently, we consider the case where mc?> ty Sn/T holds. 


Then all 63,(p) may be replaced by ordinary d(p), LE-mp) by —I’(p), and 


ae + L@--(p) by i + (m2+ 8m2(p)) +ipy a , so U[o, oo] becomes 
nee =| ds dx! A(x)  Ap(x—x!) BCx!) " 20) 
where 
FC S| Sas (dp) (3.21) 
お (ryt ) p+ C+ 5p —ho Y)/2—F€ P 


=e-2 Dt 4 p(x) 


Performing the space integral, we have 


Ula, oo] 2 (27)? ACO) BO CP —K-K)aP| é, on BOE. 
“#0 0 


( | tek +e€g/ Po)t 145 ae tat eT PMT —) ) 


; (3.22) 
Ez + €z—Po Ep + €yv—Po 


When the mass m of the intermediate particle is large, the second term is 
small enough to be neglected in comparison with the first term in this sharp 
resonance approximation. Then, 


(21 g2)? | adP dK 1 
WD eS ~ ? 
(27)5 ) 16 By Ep—ep) Ep Ep-k ed 


2 2 — 
far, (1—e 4 y+4e- 2 sin? (po ep Chan Ga EE, 3.23) 


=) Cpo—ep)? + (yC p)/2»* 
Using the formula 


(1; sin? 


BES 
a agent gon = (45 (Caos )) (3.24) 
we have 


(dP) /(27)3 
SEp(Ep—ep) en 


Wi) =2m gi (Ler) =[aP IP 1—e-"P") (3,25) 
This asymptotic time dependence of the transition probability has just the 
same form as that of (3.1). It is observed here that the effect of the detailed 
properties of the second decay is fading away exponentially from the expres- 
sion of the transition probability, and after a long time the expression tends 
to the result previously calculated from the S-matrix. Thus, we have 
obtained the divergence-free and physically acceptable result. : 


= pee Sima lala all 
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(4) several remarks 


In the same way as above, the incident line in Fig. 7 must be corrected 
as illustrated in Fig. 8. Then, if the rate of the first decay is greater than 
that of the second, the detailed nature of the first decay will fade away 
from the expression of the transition probability as if the history of the 
creation of intermediate state were completely forgotten. Such situation are 
frequently met in the theory of the nucler reactions, such as the compound 
nucleus etc., or in the case of the Hohlraumstrahlung and in the Fermi’s 
theory™ of multiple production of mesons. 

In the nuclear processes, 7r°-meson can decay into two photons through 
the virtual as well as the real state. Roughly speaking with our method, 


the decay constant ty is so small compared with the meson mass mc? that 


the decay process through virtual states will be very rare. 


Fig. 9. - nucleon line, —__~~—. meson line, ------ photon line. 
m, m', n and n! are nucleon states. A and H are their interactions. 


Next, we consider the photo-nuclear effect. It is often said that in these 
processes the mesons emitted and then re-absorbed in the nucleus play the 
important roles.) 12) At high energy, the cross-section depends mostly on 
the character of meson-production only and the detailed nature of meson- 
absorption is not important. This is just the effect of damping by absorption, 
and will be treated well according to our method. We take here only the 
simplest case. (See Fig. 9) The transformation matrix M becomes 


M = | Ap(4—x!) H(x!) = jer dt K(rt) eEn Ertl a Hir+é,0) x 


区 a ApCE, 7) eRn Em =| dr! K(r) Gr HG") (3.26) 


10) E. Fermi, Prog. Theor. Phys. 5 (1950), 570. 
11) S. Kikuchi, Phys. Rev. 85 (1951) 1062 and ibid 86 (1952) 41. 
12) R. Wilson, Phys. Rev. 86 (1952), 125. 
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where 
Gir) = whe t) e-4Et dt = —2t py le a eib,r 
も ) (273 P+ w—(4EY—il’( AE) 
| eect rez, 
~~ Aqr r (3.27) 


P=V(AEY—pe+il'(AE), [(4E)=4E-y(4E) 
Assuming the absorbing-source density is uniform in the nuclear radius R, 
and after some approximate calculation we obtain the following form). th 


vfs R 
| M |2~ = |K-H|2r~(4E) (1—e XAEB ) (3.28) 


lyway osvGd se) 


AE) 
NAE) wv AE)’ 


where aE P(AE)=V AEY— 


v4E)= 


Therefore, a picture of the free path of meson is asymptotically contained 
in this problem. 

In our above treatment we have considered only a correction of the 
internal line by the lowest order proper graph of self-energy. According to 
Dyson’s integral equation, our method consists in the replacement of ©* by 
its first term Sty in the expansion. So, this may be regarded as one-particle 
analogue of the “ladder approximation” introduced by Bethe and Salpeter! 
in their theory of the relativistic two body problem. 


+) We expand P= (AE)2—p2+il(AE) in powers of (AE) 


AE i 
tt) CAE) = AoE has some approximate meaning of (group) velocity of mesons. 


13) E.E. Salpeter and H.A. Bethe, Phys. Rev. 84 (1951), 1232... 
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